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ABSTRACT 

Uniform  and  isotropic  mathematical  models  of  the  expand- 
ing universe  usually  predict  an  initial  singularity  of  infi- 
nite mass  density  and  space  curvature.   To  study  possible 
mechanisms  which  would  avoid  the  occurance  of  these  singular- 
ities, non-uniform  cosmological  models  based  on  Einstein's 
field  equations  are  investigated  in  which  random  perturba- 
tions of  long  wave  lengths  are  superimposed  on  the  Robinson- 
Walker  metric  of  the  unperturbed  models.   Techniques  of  fluid 
turbulence  theory,  used  to  describe  random  fields  by  a  hier- 
archy of  central  moments  of  the  random  perturbations,  are 
applied  to  describe  the  dynamics  of  these  moments.   For  the 
case  of  small  perturbations  the  hierarchy  is  truncated  and 
solutions  are  found.   The  solutions  are  either  growing  or  de- 
caying perturbations  leading  to  R   extra  terms  in  the  usual 
cosmological  equations  for  the  curvature  radius  R.   The  re- 
sult agrees  with  the  small  perturbation  Fourier  series  expan- 
sion analysis  which  exists  in  the  literature.   Based  on  the 
upper  limit  of  the  anisotropy  of  the  3°  K  background  radia- 
tion, the  growing  perturbation  model  predicts  a  maximum 
expansion  even  for  k=0,  Euclidean  spaces.   The  decaying  per- 
turbation solutions  give  extra  terms  of  the  form  1/Rm  with 
m>4  in  the  cosmological  equations  and  indicate  that  the 
mechanism  of  long  wave  random  perturbations  may  prevent  the 
original  singularity  and  make  oscillatory  models  possible. 
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I.   INTRODUCTION 

Perhaps  one  of  the  most  exciting  discoveries  of  the 
Twentieth  Century  thus  far  is  Hubble 's  now  famous  discovery 
(1936)  that  the  relative  shift  in   wavelength  —-   in  the 

A 

spectral  emissions  of  each  star  is  proportional  to  its 
astronomical  distance  L  [1].   Modern  cosmology  generally 
recognizes  this  phenomenon  to  be  the  consequence  of  an 
"expanding  universe"  whose  dynamics  are  determined  by  the 
Einstein  field  equations  [2,  21].   The  concept  of  an  expand- 
ing universe  presents  various  problems  to  those  who  are 
involved  with  the  systematic  study  of  our  universe  as  an 
entity.   Because  the  term  "expanding  universe"  implies  that 
all  distances  between  all  points  are  on  the  average  increas- 
ing, it  is  natural  to  ask  if  our  universe  was  ever  in  a 
condition  in  which  all  lengths  were  zero.   It  is  also  natural 
to  ask  if  this  so-called  expansion  is  homogeneous  and  iso- 
tropic at  all  points  in  the  universe. 

If  the  answer  to  the  first  question  is  yes,  then  it  is 
apparent  that  sometime  in  the  distant  past  the  mass  density 
of  the  universe  was  infinite.   But  Wheeler  [20]  has  pointed 
out  that  such  a  state  Is  rather  dubious.   He  calls  attention 
to  the  fact  that  it  is  possible  to  construct  a  "length"  L* 
by  a  combination  of  the  gravitational  constant  G,   the 
light  velocity  c,  and  the  Planck  constant  ti  such  that 


L*  = 


■ho 


c3 


,-33 


=  1.6  x  10   cm.  (1-1) 


He  suggests  that  this  length  might  be  a  characteristic 
length  that  is  associated  with  a  quantum  theory  of  gravita- 
tion.  He  points  out  that  because  of  the  smallness  of  this 
length,  quantum  gravitational  effects,  though  negligible 
now,  could  in  principle  be  important  during  the  very  earli- 
est stages  of  the  expansion  of  our  universe  [5].   Hence, 
it  is  believed  that  the  undesirable  singularity  associated 
with  the  infinite  density  might  be  avoided.   One  reason, 
he  says,  why  such  quantum  effects  might  be  significant 
during  the  highly  contracted  phase  of  our  universe  is  that 
the  estimated  cross  section  for  pairs  of  gravitons  to  produce 
pairs  of  electrons  and  pairs  of  photons  is  comparable  when 
the  wavelength  of  the  compacted  gravitational  radiation, 
which  is  assumed  to  be  present,  is  of  the  order  of  L* .   In 
order  to  describe  such  a  situation,  it  is  clear  that  a 
quantum  theory  must  be  used  to  treat  this  obviously  non- 
classical  phenomenon  [20]. 

Various  workers  have  begun  the  task  of  incorporating 
quantum  effects  into  gravitational  theory.   Some  investiga- 
tors (Misner  [16]  and  Parker  [17]  for  example)  do  this  by 
various  canonical  procedures,  but  others  have  taken  a  differ- 
ent approach.   V/oehler  [3*0  suggests  that  it  might  be  fruit- 
ful to  postpone  a  correct  quantum  theory  of  gravitation  and 
try  a  semiclassical  approach.   A  way  he  suggests  to  do  this 
is  to  presume  that  the  field  variables  g    (called  the 


metric  tensor)  of  Einstein's  field  equations  randomly  fluc- 
tuate with  position  and  are  characterized  by  a  "mixing 
length"  L,  which  for  quantum  fluctuations  would  take  the 
value  L* .   By  considering  the  average  of  the  field  equa- 
tions, he  suggests  that  it  might  be  possible  to  avoid  singu- 
lar solutions  if  the  extra  source  term,  which  appears  as  a 
consequence  of  the  averaging  and  the  nonlinearity  of  the 
field  equations,  is  dominant  during  the  early  phase  of  high 
contraction  of  our  universe.   Ginzburg  and  other  co-workers 
[7]  use  a  similar  approach.   They  attempt  to  estimate  the 
functional  form  of  this  extra  source  term  by  expanding 
this  source  term  in  a  series  expansion  of  the  average  curva- 
ture tensor.   It  is  interesting  to  note  that  both  the  tech- 
niques of  Woehler  and  Ginzburg  are  reminiscent  of  an  earlier 
technique  used  by  Welton  [33]  to  estimate  the  Lamb  Shift 
by  a  semiclassical  model  that  also  involves  fluctuations, 
but  of  the  position  variable  of  the  electron. 

However,  the  point  in  the  above  discussion  that  is 
germane  to  the  introduction  of  this  thesis  is  that  there 
are  today  serious  considerations  of  using  a  randomly  fluctu- 
ating metric  tensor  g   as  a  prelude  to  developing  a  realis- 
tic quantum  theory  of  gravitation. 

There  is  also  another  aspect  of  cosmology  that  requires 
the  use  of  a  nonuniform  fluctuating  metric  tensor  g 
This  consideration  is  brought  out  by  the  second  question 
that  asks  if  the  expansion  is  everywhere  homogeneous  and 
isotropic.   It  is  a  matter  of  history  that  one  of  the  first 
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attempts  to  describe  the  universe  from  a  cosmological  point 
of  view  assumes  that  all  matter  and  all  space  are  homogene- 
ous and  isotropic  on  a  sufficiently  large  scale.   This 
theoretical  assertion  is  called  the  Cosmological  Principle 
and  was  formulated  in  1935  [2,  21,  22]. 

Although  the  use  of  the . Cosmological  Principle  appears 
regularly  in  the  literature,  its  criticism  is  still  very 
much  alive  today.   De  Vaucouleurs  [31] >  in  a  recent  paper, 
points  out  that  because  current  observations  have  as  yet 
failed  to  demonstrate  the  existence  of  that  sufficiently 
large  scale  beyond  which  space  is  homogeneous  and  isotropic, 
it  seems  unlikely  that  it  would  be  easily  found  (if  such  a 
length  scale  exists  at  all)  by  a  slightly  deeper  astronomical 
survey. 

Because  of  this  and  other  arguments,  many  workers  have 
been  led  to  consider  nonuniform  cosmological  models.   One 
way  this  is  done  is  to  consider  the  dynamics  of  an  anisotrop- 
ic universe  which  has  two  equivalent  and  one  inequivalent 
directions  in  space.   The  dynamics  of  such  models  can  be 
determined  exactly  by  solving  the  appropriate  field  equations 
[27].   Another  method  to  construct  nonuniform  cosmological 
models  is  to  consider  small  linear  perturbations  away  from  a 
spatially  homogeneous  and  isotropic  cosmological  metric . 
Sachs  and  Wolfe  [24]  do  this  and  are  able  to  express  these 
perturbations  as  a  function  of  time.   Their  solutions  contain 
some  arbitrary  functions  of  the  three  spatial  coordinates 
which  are  determined  by  initial  conditions  of  a  given 
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arbitrary  perturbation  field.   These  functions  then  are 
generating  functions  which  correspond  to  a  given  nonuniform 
model. 

Again  the  important  point  for  the  reader  of  this  paper 
to  note  is  that  cosmological  models  with  a  variety  of  non- 
uniformities  are  being  actively  considered  by  modern  cosmol- 
ogists.   These  nonunif ormities  may  be  zero  point  quantum 
mechanical  fluctuations  or  they  may  be  fluctuations  due 
to  a  large  scale  uneven  mass  distribution.   In  fact,  Ginz- 
burg  and  others  [7]  point   out  that  both  types  of  fluctua- 
tions are  described  by  the  same  formalism,  a  metric  tensor 
that  fluctuates  randomly  about  its  average  value.   Matzner 
[14]  points  out  that  this  formalism  is  appropriate  for  a 
universe  whose  average  background  is  uniform  but  contains 
gravitational  radiation  of  an  arbitrary  size  L.   Clearly, 
some  techniques  are  needed  which  describe  these  random  per- 
turbations of  arbitrary  magnitude  and  size  L  in  some  natural 
way. 
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II.   MOTIVATING  PHILOSOPHY  OF  THE  MODEL 

In  193*1  *  Milne  and  McCrea  chose  to  describe  our  universe 
in  terms  of  a  Newtonian  model.   This  model  treats  the  uni- 
verse as  a  streaming,  uniform  fluid.   Surprisingly,  they 
showed  that  such  a  model  is  in  many  respects  equivalent  to 
the  models  of  relativistic  cosmology  [2].   Their  model  has 
the  distinct  advantage  of  being  mathematically  much  less 
complex  than  models  developed  within  the  framework  of  gener- 
al relativity.   Another  advantage  is  that  fluid  flows  have 
been  well  studied  and  it  would  seem,  therefore,  that  some  of 
the  fruits  of  fluid  mechanics  might  be  transferable  to  the 
field  of  cosmology. 

Such  a  possibility  has,  in  part,  been  the  motivation 
for  preparing  this  dissertation.   In  particular,  the  study 
of  turbulent  fluid  flows  which  deals  with  random  perturba- 
tions of  fluid  flow  fields  would  seem,  on  the  basis  of  the 
above  discussion,  to  be  relevant  to  the  study  of  random 
perturbations  of  the  "flow  field"  associated  with  the 
universe . 

It  is  well  known  that  fluid  flows  obey  the  Navier- 
Stokes  equation.   This  equation  relates  the  velocity  field 
V  to  a  pressure  field  P  associated  with  the  flow.   It  is 
also  well  known  that  if  both  the  velocity  and  pressure 
fields  are  decomposed  into  an  average  field  plus  a  fluctua- 
tion field  such  that 
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V  =  <V>  +  v 

P  =  <P>  +  p  (2-1) 

where  <  >  denotes  "average,"  then  the  Navier-Stokes  equation 

can  be  placed  into  a  form  which  is  symbolically  identical 

-»•  ■*■ 

with  the  Navier-Stokes  equation  for  V  and  P  but  with  V 

replaced  by  <V>  and  P  replaced  by  <P>.   Owing  to  the  non- 
linearity  of  the  Navier-Stokes  equation,  additional  source 

terms  which  involve  the  correlations  <v.v.>  appear  (where 

i  J 

i  and  J  denote  the  i  and  j    components  of  the  fluctuation 
v)  [4].   These  extra  terms  are  called  the  Reynold  Stresses 
and  may  be  interpreted  as  additional  stresses  which  act 
upon  <V> . 

There  are  a  variety  of  methods  that  fluid  dynamicists 
use  to  determine  the  functional  form  of  the  Reynold  Stresses. 

An  important  way  is  to  generate  a  hierarchy  of  moment  equa- 

-> 

tions  in  v  by  multiplying  the  unaveraged  Navier-Stokes 

equation  with  all  orders  of  V  and  averaging.   Eecause  of 
the  nonlinearity ,  each  new  equation  so  generated  contains 
a  moment  of  one  order  higher  than  what  would  be  present 
if  the  Navier-Stokes  equation  were  linear.   Although  any 
number  of  moment  equations  could  in  principle  be  generated, 
there  always  exists  one  more  unknown  than  there  are  equations 
This  problem  is  called  the  closure  problem  of  turbulence 
theory.   Much  research  has  gone  into  finding  ways  to  replace 
this  indeterminate,  infinite  set  of  equations  by  a  plausible 
determinate,  finite  set  so  that  useful  information  about 
turbulant  flow  patterns  might  be  obtained  [3]-   However,  in 

1H 


spite  of  the  difficulties  associated  with  the  closure  prob- 
lem, it  should  be  pointed  out  that  the  turbulence  approach, 
which  is  a  sufficiently  general  formalism,  can  in  principle 
describe  any  sort  of  random  perturbation  field  of  any  magni- 
tude.  (A  short  summary  of  turbulence  theory  of  fluids  is 
given  in  Appendix  A.) 

The  Einstein  Field  equations  are  likewise  nonlinear 
differential  equations.   Their  field  variables  are  the 

metric  tensor  g   and  the  energy-momentum  tensor  T 

toyv  toJ  uv 

Because  certain  aspects  of  the  field  equations  of  general 

relativity  resemble  those  of  fluid  mechanics,  it  is  not 

unreasonable  to  presume  that  the  Einstein  field  equations 

could  be  expressed  in  terms  of  the  average  metric  tensor 

<g   >,  the  average  energy-momentum  tensor  <T   >,  and  some 
&yv  UV 

residual  terms  that  appear  in  a  way  analogous  to  the  Reynold 
Stresses  owing  to  nonlinearities .   It  might  also  be  possible 
to  generate  a  moment  hierarchy  which  can  be  used  for  the 
determination  of  these  extra  terms  if  an  appropriate  closure 
technique  could  be  found. 

It  is  this  speculation  therefore,  based  on  the  similar- 
ities between  fluid  mechanics  and  general  relativity,  that 
serves  as  the  motivation  for  constructing  such  a  technique 
which  might  be  useful  in  solving  perturbation  problems  of 
general  relativity.   The  remainder  of  this  dissertation  is 
a  reply  to  this  conjecture. 
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III.   INITIAL  DESCRIPTION  OF  THE  MODEL 

Given  that  it  is  possible  to  generate  a  hierarchy  of 
meaningful  moment  equations  from  the  Einstein  field  equa- 
tions, it  is  necessary  to  truncate  them  in  some  manner. 
Then  a  solution  for  the  "Reynold  Stresses"  can  be  determined 
which  will  be  found  to  contain  second  order  moments  in  the 
perturbed  metric  and  its  derivatives.   Just  as  in  fluid 
turbulence  theory,  the  closure  of  the  moment  hierarchy  is 
nontrivial  because  it  requires  special  restrictive  assump- 
tions concerning  the  nature  of  the  random  perturbation 
field  and  does  not  just  "fall  out"  of  the  theoretical  fra 
framework. 

In  order  to  see  how  this  truncation  might  be  accomplished, 
consider  first  the  following  three  important  truncation 
procedures  which  are  used  in  fluid  turbulence  theory  to 
close  its  unlimited  set  of  equations  which  contain  moments 
of  all  orders  of  the  velocity  v  [3]. 

1.  Phenomenological  approximations  -  These  involve  postu- 
lating (on  dimensional  grounds  or  by  intuition)  an 
extra  relation  between  various  moments  to  bring  the 
number  of  equations  up  to  the  number  of  unknowns. 

2.  Higher-order-moment-discard  approximation  -  This  tech- 
nique simply  discards  the  higher  order  moments  because 
of  smallness  or  some  other  reason.   This  reduces  the 
number  of  unknowns  to  the  number  of  available  equations. 
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3.  Cumulant-discard  approximation  -  This  procedure  relates 
a  higher  order  moment  to  products  of  lower  order  moments 
and  thereby  effectively  establishes  an  extra  equation. 

All  these  closure  methods  are  used  in  fluid  turbulence 
theory  and  so  all  might  conceivably  be  acceptable  candidates 
for  closing  a  moment  hierarchy  in  general  relativity  as  well. 

If  it  is  desired  to  describe  cosmological  models  with 
a  "small"  degree  of  anisotropy  and  inhomogeneity ,  then  it 
would  seem  that  the  appropriate  truncation  would  be  number 
two,  the  higher-order-moment-discard  approximation.   This 
choice  is  reasonable  because  it  would  be  expected  that  a 
moment  of  order  n+1  would  be  negligible  compared  to  a  moment 
of  order  n.   Such  weakly  perturbed  models  are  considered, 
for  example,  by  Sachs  and   Wolfe  [24]  as  previously  stated. 

However,  not  all  cosmological  models  are  adequately 
described  by  small  perturbations.   In  a  recent  paper  by 
Woehler  [34],  a  model  is  presented  which  assumes  quantum 
fluctuations  of  the  metric  tensor  to  be  as  large  as  the 
average  metric  tensor  itself.   Certainly  the  higher-order- 
moment-discard  approximation  is  not  adequate  to  apply  to  this 
model  because  all  moments  of  all  orders  are  important,  espe- 
cially near  the  singularity.   For  this  type  of  model,  the 
first  or  third  truncation  alternative  might  be  useful. 

Therefore,  an  attempt  to  construct  a  general  perturba- 
tion technique  based  on  a  moment  hierarchy  scheme  which  is 
able  to  treat  fluctuations  of  all  shapes  and  sizes  is  dis- 
couraged by  the  fact  that  the  method  of  truncation  depends 
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upon  the  type  of  fluctuation  that  is  involved.   In  fact, 
even  the  solutions  of  the  moment  hierarchy  could,  in  general, 
be  quite  different  owing  to  the  fact  that  different  moment 
equations  could  result  from  different  methods  of  truncation. 

So  it  seems  that  a  reasonable  way  to  proceed  is  to 
consider  a  model  with  a  certain  "kind"  of  perturbation 
field  and  then  form  and  truncate  its  moment  hierarchy  by  a 
technique  that  is  compatible  with  the  type  of  fluctuations 
assumed.   Thus,  a  two  fold  purpose  would  be  achieved.   First, 
a  model  would  be  constructed  which  describes  a  special  non- 
uniformity  about  which  conclusions  can  be  drawn  that  may  be 
of  value  to  cosmology.   Second,  and  perhaps  more  signifi- 
cant, would  be  the  demonstration  that  this  procedure  does 
work  and  is  therefore  relevant  to  perturbation  problems 
of  general  relativity. 

In  order  to  best  serve  the  process  of  demonstrating 
the  relevancy  of  the  moment  hierarchy  technique  to  general 
relativity,  a  simple  nonuniform  model  is  used.   It  is  simple 
in  the  sense  that  only  one  function  is  required  to  describe 
the  perturbation  field  rather  than  ten  which  are  needed  to 
treat  a  general  perturbation  problem.   This  random  function 
is  also  assumed  to  contain  random  perturbations  which  are 
small  compared  to  the  average  field.   Thus,  the  previously 
discussed  higher-order-moment-discard  approximation  will 
be  used  for  the  truncation  of  the  moment  hierarchy  which 
will  be  constructed  in  a  way  analogous  to  fluid  turbulence 
theory.   Even  for  this  relatively  simple  model,  the 
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mathematical  manipulations  are  quite  extensive.   It  will  be 
shown  that  this  particular  model  is  relevant  to  the  study 
of  long  wavelength  perturbations  and  deformations  of  the 
size  which  approaches  the  length  characteristic  of  the 
universe  itself.   Thus,  the  model,  though  highly  restrictive, 
has  a  special  relevancy  and  direct  applicability  to 
cosmology. 
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IV.   METRIC  TENSOR 

A.   SPECIAL  FORM  FOR  THE  MODEL 

The  process  of  building  the  model  begins  with  the  formu- 
lation of  the  metric  tensor  g   .   This  can  be  accomplished 
by  giving  first  the  equation  for  the  square  of  the  line 
element  ds  relative  to  a  given  coordinate  system  x   and  then 
computing  g   (x  )  by  the  relation 

ds2  =  g  dxydxv  .  (4-1) 

To  this  end,  consider  first  the  line  element  for  uniform 
cosmological  models .   It  is  called  the  Robertson-Walker  line 
element  and  in  spherical  coordinates  is  given  by  [1] 

ds2  =  (dx°)   -  eG(x°'r)da2  (4-2) 

where 

da2    =   dr2    +   r2d62    +   r2sin20d<{>2  (4-3) 


and 


2 
.2      "I 


2 

r  „ 


k  r' 
*  ~7 


ro 


The  parameter  r0  is  an  arbitrary  normalization  constant 
with  dimensions  of  length  for  the  radial  coordinate  r. 
k  is  the  curvature  characteristic  which  may  take  only  the 
values  -1,  0,  and  +1.   The  time  coordinate  is  x   and  the 
three  spatial  coordinates  are  r,  e,  and  $ .   The  function 
R_(x°)  is  called  the  scale  length. 

u 
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Adler  [1]  points  out  that  the  spatial  part  of  the  line 
element  of  Eq.  (4-2)  describes  a  three  dimensional  geometry 
which  is  homogeneous  and  isotropic  with  a  uniform  radius  of 
curvature  RQ .   Associated  with  the  three  values  of  the 
curvature  characteristic  k  are  three  intrinsically  differ- 
ent kinds  of  uniform  geometries.   The  first,  for  which  k 
is  equal  to  zero,  is  the  familiar  Euclidean  geometry.   This 
geometry  is  based  upon  the  postulate  that  through  any   point 
outside  of  a  straight  line,  one,  and  only  one  parallel  line 
can  be  drawn.   The  case  where  k  is  equal  to  +1  corresponds 
to  the  uniform  geometry  first  discussed  by  Riemann  which 
is  based  on  the  postulate  that  any  two  straight  lines  in  a 
plane  always  intersect  each  other.   Finally,  the  case  where 
k  is  equal  to  -1  corresponds  to  the  noneuclidean  geometry 
of   Bolyai-Lobachevaki .   This  uniform  geometry  is  based 
upon  the  postulate  that  through  a  point  outside  of  a  straight 
line,  an  unlimited  number  of  straight  lines  can  be  drawn 
which  do  not  intersect  the  given  line  and  which  lie  between 
two  straight  lines  that  do  intersect  the  given  line  at 
infinity  [26]. 

Now,  if  the  function  G(x°,r)  were  replaced  by  an  arbi- 
trary function  M(x°  ,r  ,0  ,<j> )  ,  the  resulting  line  element  would 
no  longer  be  necessarily  uniform.   In  fact,  a  positive 
function  RM(x° ,r , 0 ,$ )  could  be  defined  such  that 

M(x°,r,0,4>)  _   RM  (x0>r>9>4>> 


r  2 

o 


1  + 


(4-5) 

k  r2  ^2 
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Thus,  at  every  point  xy,  a  tangent  Robertson-Walker  geometry 
could  be  constructed  with  a  radius  of  curvature  R_  equal  to 

u 

the  function  R..(x  ) .   Hence,  the  function  R„  could  be  viewed 
M  M 

in  this  sense  as  a  "local  radius  of  curvature"  or  "local 
scale  length." 

As  stated  previously,  one  of  the  objectives  of  this 
dissertation  is  to  demonstrate  the  applicability  of  a  prin- 
ciple of  fluid  turbulence  to  perturbation  problems  of  general 
relativity  by  means  of  a  simple  model.   Clearly,  a  model 
whose  geometry  is  described  by  the  line  element  of  Eq .  (4-2) 
with  G  replaced  by  M  is  simple  because  the  entire  nonuni- 
formity  is  described  by  the  single  function  M  (in  general, 
ten  independent  functions  are  needed  to  describe  a  nonuni- 
form metric  guv(x  )  because  g    is  symmetric  with  respect 
to  u,v).   Thus,  such  a  model  is  a  convenient  choice  for 
such  a  demonstration.   It  is  therefore  assumed  that  the 
line  element  of  the  model  is  given  by 

2        lUT/  „  o 


ds 


»  =  (dx°)   -  eM(x  .*.M>da.  .  (4_6) 


Because  this  model  is  to  be  described  by  statistical 
methods,  the  value  of  the  function  M  at  any  given  point 
must  possess  a  certain  amount  of  randomness.   That  is,  it 
must  be  viewed  as  a  stochastic  process  or  random  field. 
Hence,  M  is  a  function  not  only  of  the  four  space-time 
coordinates  (x°,r,6,<t>)  but  also  of  an  ensemble  parameter  K 
which  labels  each  of  those  functions  M  that  are  possible 
for  a  given  experiment.   Hence,  Eq .  (4-6)  describes  not  just 
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one  geometry,  but  an  entire  family  of  geometries  each  of 
which  must  obey  the  Einstein  field  equations. 

The  notion  of  an  experiment  in  the  study  of  stochastic 
processes  is  well  defined.   See,  for  example,  Papoulis  [18]. 
However,  for  the  purposes  here,  there  is  an  interesting 
connection  between  an  experiment  (as  used  in  statistical 
theory)  and  a  modern  concept  called  superspace.   Wheeler 
[5]  defines  superspace  as  the  manifold  each  of  whose  "points" 
represents  one  three-geometry.   A  given  submanifold  of 
superspace  is  therefore  a  history  of  the  geometry  undergone 
by  physical  space.   Hence,  the  term  experiment  as  used  in 
this  paper  would  indicate  an  observation  of  the  "state 
vector"  of  a  possible  three-geometry  in  superspace  at  a 
given  time.   Certainly,  the  set  of  possible  state  vectors 
for  a  given  experiment  is  severely  limited  by  Eq .  (4-6). 

If  the  tensor  form  of  the  line  element  given  by  Eq. 
(4-1)  is  compared  with  the  line  element  given  by  Eq .  (4-6) 
with 

x°  =  x° 

x1  =  r 

x2  =  e 

x3  -  4  ,  (4-7) 

it  is  easily  found  that  the  nonvanishing  elements  of  the 
metric  tensor  for  the  model  must  be 

600  = 

g„  =  -  eM  (4-8) 

cont . 
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2     M 

g22     =   -   re 

g33     =   -  r^sin26e        .  (4-8) 

B.   DECOMPOSITION  OF  THE  METRIC  TENSOR  INTO  ITS  AVERAGE 

AND  FLUCTUATING  COMPONENTS. 

In  fluid  turbulence  theory,  each  of  the  random  field 
variables  (pressure  and  velocity)  of  the  governing  Navier- 
Stokes  equation  is  decomposed  into  its  average  and  fluctuat- 
ing components  as  described  in  Chapter  II.   This  same 
decomposition  may  be  performed  upon  the  random  metric  tensor 
g   which  is  the  field  variable  of  general  relativity. 
That  is,  g   may  be  written  as 

g   =  <g  >  +  f  (4-9) 

&uv    ^uv     uv 

where  <g   >  is  the  average  of  g   and  f   is  its  fluctuating 
touv  °     &uv      yv 

component.    This  average  is  understood  to  be  an  ensemble 
average.   That  is,  to  take  an  average  at  each  point  x^1 , 
the  value   of  the  random  field  from  each  member  of  the 
ensemble  of  possible  states  of  the  system  which  occurs 
with  a  certain  probability  are  summed  and  divided  by  the 
number  of  ensemble  members.   As  this  number  tends  towards 
infinity,  this  sum  tends  towards  the  average  value.   (The 
concept  of  average  is  discussed  in  more  detail  in  Appendix  B. ) 

As  seen  in  Eq .  (4-8),  the  randomness  of  g   is  contained 

M 
in  the  stochastic  process  e  .   This  process  can  also  be 

decomposed  at  each  point  in  space-time  into  two  functions 

such  that 
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M     M 
e   =  <e  >  +  fM  (4-10) 

M  M 

where  <e  >  is  the  average  of  e   and  f   is  the  fluctuation  of 

M 

e   about  its  average.   By  Eq.  (4-8),  the  nonvanishing  com- 
ponents of  <g   >  and  f   must  be 
yv       yv 

<goo  >  =  1 

M 
<gu  >  =  -  <e  > 

2   M 
<g22  >  =  -  r^<e  > 

9   .   9      M 

<g33  >  =  -  r  sini<e  > 


f„  =  -  f 


M 


f  22  =  "  r2fM 


fsa  =  -  r2sin20fM   .  (4-11) 

Another  stochastic  process  t(x^)  can  be  defined  as  the 

M 

ratio  of  the  fluctuation  of  e   to  its  average.   That  is, 

T   S  -$-   .  (4-12) 

<e  > 

If  f    is  expressed  in  terms  of  t.  its  nonvanishing  elements 
yv 

become 

«         M 
fn  =  -  <e  >t 


f„  =  -  r2<e  >x 


f33  =  -  r2sin2e<e  >i  .  (4-13) 

Because  the  average  of  f..  must  be  zero  according  to  Eq.  (4-10), 
it  follows  that 

<T>  =  0   .  (4-14) 
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Also  by  Eq.  (4-10),  it  follows  that 

t  >  -  1  (4-15) 

M  M 

because  e   and  therefore  <e  >  must  be  positive.   Certainly 

models  with  small  perturbations  would  be  ones  for  which  the 

ratio  of  the  nonvanishing  components  of  f    to  <g   >  is 

yv     &yv 

small;  that  is,  for  these  models,  the  magnitude  of  x  would 
be  very  small  compared  to  one.   As  stated  in  Chapter  III, 
it  is  this  kind  of  perturbation  that  will  be  considered  in 
this  work. 

1.   Requirement  for  the  Average  Component  of  the  Metric 
Tensor 

Suppose  it  is  assumed  that  the  average  metric  tensor 
<g   >  is  Robertson-Walker.   Then,  the  average  geometry  would 
be  homogeneous  and  isotropic  everywhere .   This  assumption, 
which  is  used  extensively  in  the  literature  [7,  14,  24], 
is  convenient  because  of  its  relatively  simple  mathematical 
form.   Furthermore,  in  the  limit  of  zero  fluctuations,  the 
model  becomes  determinate  and  the  geometry  is  automatically 
uniform. 

Thus,  the  nonvanishing  components  of  the  metric  tensor 
<g   >  for  the  model  are  given  by 

<goo  >  =  1 

G(x°,r) 
<gn  >  =  -  e     » 

2  G(x° ,r) 
<g22  >  =  -  r  e     » 

<g33  >    =   _   r2sin20eG(x°>r)       .  (4-16) 
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Because  the  average  metric  is  assumed  to  be  everywhere  homo- 
geneous and  isotropic  and  is  given  by  the  above  equations, 

M      C 
it  follows  that  <e  >  =  e  .   By  virtue  of  Eqs.  (4-4)  and 

(4-5),  it  can  be  easily  seen  that 


2 


<RMZ  (x°  ,r,6,4,)>  =  RQ  (x°)  (4-17) 

where  RM(x°  ,r  ,e  ,<J>)  is  the  already  discussed  local  scale 
length  of  the  random  metric  and  R„  is  the  universal  scale 
length  of  the  average  metric.   If  AR2  (x  °  ,r ,  6  ,<J> )  is  the 

2  2 

fluctuation  of  RM   from  R„   such  that 

RM2  =  RG2  +  AR2  ,  (4-18) 


then  by  Eq .  (4-12),  it  follows  that 


AR2 

2 

RG 


(4-19) 


Thus,  t  can  be  interpreted  as  the  ratio  of  the  fluctuation 
of  the  local  radius  of  curvature  squared  to  the  average 
radius  of  curvature  squared. 

2 .   Requirement  for  the  Fluctuating  Component  of  the 

Metric  Tensor 

In  fluid  turbulence  theory,  a  convenient  require- 
ment concerning  the  statistics  of  the  perturbation  field  is 
the  requirement  of  stationarity  [13]-   If  this  statistical 
quality  is  borrowed  and  applied  to  the  present  model,  some 
simplifications  can  be  expected  by  assuming  the  statistics 
of  the  turbulent  field  t  to  be  stationary  with  respect  to 
a  change  of  coordinates  on  a  spatial  hypersurface  of  const 
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x°.   Clearly,  the  statistics  may  not  be  assumed  to  be 
stationary  with  respect  to  the  time  coordinate  x°  because 
this  would  imply  a  model  universe  that  is  not  in  a  state 
of  dynamic  evolution  (i.e.  expanding).   In  Appendix  C,  the 
implication  of  this  concept  is  discussed  and  the  important 
formula  given  below  is  derived. 

9.<h[T(xy)]>  =  0  (4-20) 

J 

where  j  =1,  2,  3  and  h(i)  is  any  function  of  the  stochastic 
process  T.   Clearly  this  is  a  very  useful  formula  for  by  it 
many  bothersome  terms  are  eliminated. 
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V.   EINSTEIN  FIELD  EQUATIONS  FOR  THE  FLUCTUATING  METRIC 

A.   DESCRIPTION  OF  THE  FIELD  EQUATIONS  IN  GENERAL 

Just  as  the  Navier-Stokes  equation  relates  the  velocity 

field  to  the  pressure  field  in  a  fluid,  so  do  the  Einstein 

field  equations  relate  the  geometry  of  space  to  the  energy 

content  of  that  space  [1].   Being  a  tensor  theory,  the 

field  equations  are  expressible  in  a  form  independent  of  any 

coordinate  system.   The  geometrical  variable  of  the  field 

equations  has  already  been  identified  as  the  metric  tensor 

g   .   The  energy  content  variable  is  called  the  energy 

momentum  tensor  T   which  is  developed  and  discussed  later. 

yv 

It  should  be  stated,  however,  that  the  field  equations  for 

the  model  itself  are  not  expressed  in  covariant  form  because 

a  special  coordinate  system  given  by  Eq .  (^1-7)  has  already 

been  imposed. 

Starting  with  the  metric  tensor  g   ,  the  field  equations 
°  yv 

can  be  constructed  as  follows.   (This  procedure  is  given  by 

Adler  [1]  and  is  used  in  the  subsequent  sections  to  build 

the  field  equations  specific  to  the  metric  tensor  of  the 

proposed  model.) 

First,  the  metric  tensor  g   is  substituted  into  the 
3  yv 

following  expression. 

[yv,3l  =  \   (3vgp3  +  9ygv3  "  Vyv)  (5-D 


V    ^ 


where  3        .   [yv,$]  is  the  Christoffel  symbol  of  the 

29 


first  kind  and  is  a  set  of  6H    elements  which  are  symmetric 
with  respect  to  their  y,  v  indices. 

Next,  a  quantity  r"   ,  the  Christoffel  symbol  of  the 
second  kind,  which  also  is  a  set  of  6k   elements,  is  calcu- 
lated by 

ryv  =  Sa3[yv,6]  (5-2) 

where  the  identical  indices  indicate  a  summation  over  all 

Q, 

3.  T  is  also  symmetric  with  respect  to  its  lower  indices. 
The  object  g  ,  called  the  second  rank  contravariant  metric 
tensor,  is  related  to  the  metric  tensor  g   by 

g  Ey3  ■  6y   .  (5-3) 

Next,  the  Ricci  tensor  R   is  formed  by 

r   =  8  ra  -  a  rp  +  ra  rp  -  ra  rp   .     (5-*0 

yv    v  pa     p  yv     pv  ya     pa  yv 

It  is  a  set  of  16  elements.   The  Ricci  tensor  is  then  used 
to  calculate  the  mixed  Ricci  tensor  R   given  by 

KXv  =  SX\v  •  .  (5-5) 

Next,  the  Ricci  scalar  R  is  formed  by  the  contraction 
of  the  mixed  Ricci  tensor  R   given  by 

R  =  RXX  (5-6) 

Finally,  the  mixed  Ricci  tensor  of   Eq.  (5-5)  and  the 
Ricci  scalar  of  Eq.  (5-6)  are  used  to  form  the  mixed  Einstein 

tensor  G 

v 
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G*   =  RX   -  \  6X    R   .  (5-7) 

The  Einstein  field  equations  are  then  formed  by  equating 

the  mixed  Einstein  tensor  G    to  another  tensor  kTX   called 

v  v 

the  mixed  energy  momentum  tensor,  an  object  which  is  related 

to  T   and  is  discussed  later, 
yv 

Thus,  the  mixed  form  of  the  Einstein  field  equations 
as  they  are  used  in  this  work  are 


where 


kTAv  =  GXv  (5-8) 


~  2 


B.   CONSTRUCTION  OP  THE  CHRISTOFFEL  SYMBOLS  OF  THE  SECOND 

KIND  FOR  THE  MODEL 

The  task  of  constructing  the  field  equations  for  the 
proposed  model  itself  may  now  begin  as  outlined  in  the 
preceeding  section.   From  the  field  equations,  the  moment 
hierarchy  will  be  determined. 

An  important  object  that  is  immediately  needed  to  form 

the  Christoffel  symbols  is  the  second  rank  contravariant 

metric  tensor  g   .   This  object  can  be  computed  from  Eqs . 

(4-8)  and  (5-3).   The  nonvanishing  elements  of  this  tensor 

are 

g00  =  1 

n       -M 
g   =  -  e 

=  _e^  (5_9) 

s       2  cont. 

r 
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-e 


-M  cont. 

g33  =  r (5_9) 

r2  sin2e 

The  next  step  is  to  substitute  the  metric  tensor  of  Eq .  (4-8) 

into  Eq.  (5-1)  which  yields  [uv,3],   the  Christoffel  symbol 

of  the  first  kind.   Because  the  second  rank  contravariant 

metric  tensor  g^   is  diagonal,  the  calculation  of  V        ,  the 

Christoffel  symbol  of  the  second  kind,  by  means  of  Eq.  (5-2) 

is  simplified. 

The  calculation  of  r    is  still  somewhat  lengthy  but  is 

yv 


entirely  straightforward.   The  nonvanishing  components  of 
yv 


Ta   turn  out  to  be 


o     1   M 
rn   =  |  eu30M 

r^   =  ^r2eM9oM 


0      1  M 

r33      =   |  r2sin2ee    80M 
foi      =    fio      =    2    ^o^ 


I'll      =    i    9   M 

11  2      r 


r*     =  -  |r2OrM  +  I) 

T*      -   -   i  r2sin2eO    M  +   |) 


2  2  1 


r0     -  r»     -  i  90M  (5-io) 


1 02  "     *  20  ~     2  ° 


cont . 
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2  ii  cont. 

Tu       "    "    2   T  9GM  (5-10) 


r2 


fa      =    r j      =   \   (3pM   +  |-) 


r^     =  \ 


2 


2  2  1.. 

Is-       -   Ta      -   |   a^M 

T33      =    -   |   sin26(2    cote    +    3   M) 


3  3  1 

r03     =  r^     =  ±  30m 


ru3     =  -  I  — i a  m 

^   r2sln2e      * 


r3?     =  r133     ■  |  (3rM  +  |) 


3  l      l 

2    sin2e      * 


r*2      =    Ta      s|   (2    cote    +    3.M) 


2 


3  1 

2    <V 


r33     =  -  8,m 


C.   CONSTRUCTION  OF  THE  RICCI  TENSOR  FOR  THE  MODEL 

The  next  step  is  the  calculation  of  the  Ricci  tensor 
R   from  Eqs.  (5-*0  and  (5-10).   However,  it  is  anticipated 

yj,   v 

that  only  the  diagonal  elements  of  the  Ricci  tensor  are 
needed.   (This  will  be  justified  in  the  subsequent  section.) 
These  diagonal  elements  are 

Roo  ■  |  302M  +  }   O0M)2 

_  a  1    1         2  i  1  2  (5-H) 

Ru     =    3pM   +    2   —   9QM  +    2      2     .    2a    V  cont. 

r  r    sini      Y 
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-ie»3^M-|eM(30H)2  +  |i-OeM) 


r2 


r  r  sin  0    y 


1        M  1        2         2 

R22  =  -  ^  r2e  3jM  +  =r   r23  M  +  8nM 
2       0     2     r      0 


+  \   ■-—-  9IM  +  F  r'C^M)2  +  ^   — 


.n2 


sin  9 

+   I  cot©    90M   +   I  r3   M   -   |  r2eM(30M)2 


R33     =   -   I  r2sin20eM3*M  +   |  r2sin2692M   +   |  sin2090M 
+    3|M  -  I  r2sin20eMOoM)2    +    sin2©    cot0    dQM 
+   ^  r2sin20(8rH)2    +   J-  sin20(3eM) 


+  I  r    sin208    M    .  (5-11) 


D.   CONSTRUCTION  OF  THE  DIAGONAL  ELEMENTS  OP  THE  MIXED 

EINSTEIN  TENSOR  FOR  THE  MODEL 

In  order  to  construct  the  mixed  Einstein  tensor,  the 

two  objects  R   and  R  must  be  found, 
v 

The  mixed  Ricci  tensor  R   is  easily  calculated  from 

v 

Eq.  (5-5).   Since  g    is  a  diagonal  tensor,  R   can  be 
written  as 


Rx  =  (x)(x)R  (5_l:L) 

v  kX)v 

where  the  parenthetical  indices  indicate  no   summation  with 
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respect  to  A.   Equation  (5-11)  justifies  the  assertion  that 

only  the  diagonal  components  of  R   are  needed  if  only  the 

yv 

diagonal  elements  of  R   are  to  be  calculated  because  no 

v 

off -diagonal  mixing  occurs  by  virtue  of  the  diagonality 
of  g 

The  Ricci  scalar  R  is  easily  calculated  from  the  diago- 
nal elements  of  R   by  Eq .  (5-6).   Finally,  the  components 

of  the  mixed  Einstein  tensor  G   given  by  Eq.  (5-7)  turn 

v 

out    to   be 

G°o    =    -   }   (30M)2    +   e~M[^M   +   i-   3*M 

+ 3*M   +   i  —   (3nM)2    +   \   3    M 

r2sin2e      *  *   r2        °  r     r 

+   S^li   3flM   +   i O.M)2    +   i    (3    M)2] 

r2         6  *   r2sin26        *  ^        r 

G1!    -   -    ajM  -  f   O0M)2    +   e~M[   ^   3*M 

+  1  1 32M  +  i  3   M  +  i  £2£!  3AM 

2   r2sin2G      *  r      r  2      r2         6 

+  If   OrM)2] 
0«t    =   -3jM-3    (3oM)2    +  e-ML   1   a.M  +   1  _i_  3JM 

r  r 


G33    .    _    3^M_    3    (3oM)2     +   e-Hc    1    a,M+   |i-32M 


(5-12) 
cont . 
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11  in  2  cont. 

^F^UT^C,")        •  (5-12) 

r2  sin2  6        T 

The  off-diagonal  components  are  not  given  because  they 

are  not  needed.   However,  calculations  similar  to  the  ones 

required  to  obtain  Eq .  (5-12)  reveal  that  the  off -diagonal 

components  of  GA   are  in  general  non  zero  and  contain  mixed 

v 

partial  derivatives  of  second  order  and  second  degree  prod- 
ucts of  first  partial  derivatives  with  respect  to  all  four 
coordinates . 

It  should  be  noted  that  if  M  is  a  special  function  of 
only  the  coordinates  x°  and  r,  namely 

M  =  G(x°,r)  (5-13) 

where  G  is  given  by  Eq.  (4-4),  then  the  line  element  is 
identical  to  the  line  element  of  uniform  cosmological  models 
[1].   In  this  case,  as  should  be  expected,  Eq.  (5-12)  is  in 
fact  identical  to  the  diagonal  elements  of  G   for  these 
cosmologies . 

E.   CONSTRUCTION  OP  THE  DYNAMIC  EQUATIONS  FOR  THE  MODEL  VIA 

THE  EINSTEIN  TENSOR 

It  should  be  noted  that  only  two  "dynamic  equations" 
for  g   are  used  to  describe  uniform  cosmological  models[l]. 
But  each  expression  of  Eq.  (5-12)  represents  one  side  of 
four  "dynamic  equations."   If  the  proposed  model  of  this 
study  is  to  pass  to  a  uniform  cosmological  model  for  arbi- 
trarily small  perturbations,  then  the  number  of  equations 
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which  should  be  treated  as  dynamical  equations  ought  to  be 
two  rather  than  four. 

In  order  to  see  how  to  construct  two  equations  from 
equation  set  (5-12),  it  should  be  observed  that  G  0  contains 
second  partial  derivatives  with  respect  to  all  four  coordi- 
nates x°,r,0,  and  <j> .   But  each  of  the  remaining  objects  G1  , 
G22,  and  G33  contain  second  partial  derivatives  with  respect 
to  only  two  of  the  three  possible  spatial  coordinates. 
This  assymetry  can  be  removed  if  G°   is  taken  as  one  side  of 
one  of  the  desired  two  'dynamical  equations"   and  one  third 

the  sum  G^  .   such  that 
J 

GJ  .  =  G1   +  G2,  +  G3  (5-1*0 

J        1        2        3 

is  taken  as  the  other.   These  two  equations  contain  dynam- 
ical information  concerning  g   and  contain  second  partial 

yv 

derivatives  with  respect  to  all  coordinates.   They  are  also 
found  to  degenerate  into  the  usual  two  cosmological  equa- 
tions when  G  is  replaced  by  M  [1]. 

By  virtue  of  Eq.  (5-8),  the  two  dynamic  equations  for 
the  model  are 

kT°0  -  G°0  (5-15) 


and 

E   t-3'  .  =  i-  (5-16) 


GJ 


From  Eq.  (5-12),  G^ ,  is  given  by 

j 

GJ   =  _  3a2M   9  (g  m)2  +  e~M[9*M  +  i-  din  (5-17) 

J        o     4    o  r    r2        cont> 
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1  -i  2-1-1  2      cont. 

+  0?M)    +   £   0    M)      +   ^  ±-   (3    M)         (5-17) 

r2  sin20        ^  M  q   r2        u 

+  £ O.M)2    +  §   (3    M)    +  22M.  3    M]. 

4   r2sin2e        *  r        r  r2         6 

For  convenience,  G°   can  be  expressed  in  terms  of  GJ' .  . 

0  J 

G\    =    G3.    +  382M  +  |  (8oM)2  .  (5-18) 


F.   CONSTRUCTION  OF  THE  ENERGY-MOMENTUM  TENSOR  FOR  THE  MODEL 

To  proceed  further,  it  is  necessary  to  express  the 

energy-momentum  tensor  in  a  form  that  is  specific  to  the 

metric  used  for  this  model.   It  is  to  be  understood  that 

the  mixed  energy-momentum  tensor  T^   is  also  a  stochastic 

v 

process' because  it  depends,  via  the  Einstein  field  equations, 

upon  the  metric  tensor  which  is  a  stochastic  process.   The 

general  form  of  the  energy-momentum  tensor  is  known  to 

depend  upon  three  quantities,  D  the  proper  mass  density, 

P  the  proper  pressure,  and  u  the  velocity  of  the  medium. 

These  quantities  D  and  P  are  termed  "proper"  because  they 

are  referred  to  the  proper  frame  which  is  that  frame  that  is 

instantaneously  moving  with  the  medium.   From  any  other 

frame,  the  well  known  effects  upon  the  mass  and  volume  due 

to  the  relative  velocity  u  must  be  taken  into  account. 

The  representation  of  the  general  energy-momentum  tensor 

with  respect  to  the  metric  that  is  used  in  this  model  is 

derived  in  Appendix  D.   All  of  the  elements  of  this  tensor 

are  given  by  Eq.  (D-19).   The  two  related  objects  which  are 

important  here  are  T°   and  TJ .  ; 

o       J 
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P     2 

To   =  d  +  —  u 
c2 


iJ    -     n„2     P 


where 


TJ   =  -  Du2  -  —  (3  -  2u2)  (5-19) 

J  c2 


D  =  proper  density  D  x  - 


1-u2 

P  =  proper  pressure  P  x  (5-20) 

1-u2 

u  =  local  velocity  of  the  matter  in  units  of  c. 

Because  of  the  recent  discovery  of  a  possible  3°K  black- 
body  background  radiation  [6,  23],  a  relevant  energy-momentum 
tensor  is  one  which  describes  disordered  electromagnetic 
radiation  of  this  type.   However,  this  radiation  is  believed 
to  be  no  longer  important  to  the  dynamics  of  the  universe 
(because  now  the  matter  density  far  exceeds  the  radiation 
density)  although  it  is  believed  to  have  been  in  the  more 
distant  past.   The  energy-momentum  tensor  which  is  relevant 
to  the  current  stage  of  expansion  is  one  which  describes  the 
random  motion  of  matter  in  the  form  of  galaxies  and  clusters 
of  galaxies.   Eq.  (D-1))  of  Appendix  D  gives  for  these  two 
cases 

P  =  ac2D  (5-21) 

where   a  =  0  if  universe  is  matter  dominated  and 

a  =  1/3  if  universe  is  radiation  dominated. 

Hence,  by  Eq.  (5-20) , 

P  =  ac2D.  (5-22) 
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VI.   MOMENT  EQUATION  HIERARCHY 

A.   AVERAGE  OF  THE  DYNAMIC  EQUATIONS 

Equations  (5-15)  and  (5-16)  are  the  "equations  of  motion" 
for  the  assumed  metric  tensor  of  Eq.  (4-8)  and  the  energy- 
momentum  tensor  appropriate  for  that  metric.   These  equations 
should  be  viewed  as  the  analogue  of  the  Navier-Stokes  equa- 
tion  for  V  and  P  given  by   Eq .  (A-l)  because,  as  yet,  the 
decomposed  metric  tensor  of  Eq .  (4-9)  has  not  been  inserted 
into  them.  s 

However,  the  substitution  of  the  decomposed  metric 
can  be  accomplished  rather  easily.   A  quick  inspection  of 
GJ .  and  G   reveals  that  they  are  made  up  exclusively  of 

partial  derivatives  of  the  stochastic  process  M  and  a  factor 

-M 
e    .  From  Eqs.  (4-10)  and  (4-12),  it  is  easily  found  that 

eM  =  eG(l  +  t)  (6-1) 

where  e   is  associated  with  the  average  metric  and  x  with 
the  fluctuating  part  of  the  metric. 

Equation  (6-1)  can  be  used  to  separate  GJ .  and  G  0  into 
two  parts,  one  which  contains  only  the  nonrandom  function 
G  and  the  other  which  contains  the  random  function  x  in 
combinations  of  G  and  x.   The  part  which  is  made  up  only  of 


G  terms  must  be  G  .  and  G  „  based  on  the  average  metric 

j 

(because  this  part  cannot  depend  on  x)  while  the  remaining 


part  must  be  the  source  terms  (because  these  terms  must 

o        i 
depend  on  x) .   The  resulting  expressions  for  G   and  G  . 

j 
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are  (using  Eq.  (4-4)  which  is  the  definition  of  eG  and  with 
R  equal  to  RQ  and  Eqs .  (5-15)  and  (5-16)) 


0        0 


-G 


3j£  _  3R1 
R2 


1+T 


1  +  T 


R- 


3  (V> 


3R' 
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1  +  T 
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2  I 
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r 2sin2 


2  1 
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3  ^ill 

7T  — _ 

4  (1+t) 


cote 


3i£ 
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1  + 


kr2 

4r2 
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V 


1+T 


1  +  T 


3  <V> 


4  2 

4   (1+T)' 


(6-2) 


kTj\  =  Gj.  =  g°  -  ^~ 

J      J      o     R 


33Jt   3  O0t) 

1+T+^(1+T)2 


6R_ 
R 


V 


1+T 


(6-3) 


If  these  equations  are  averaged,  the  result  is  the 

desired  dynamic  equations  for  the  average  variables  <g   > 

and  <T°  >. 
o 

These  averaged  equations  can  be  immediately  simplified 

by  using  the  assumption  of  spatial  stationarity  that  is 

discussed  in  Ch.  IV.   Thus,  all  averages  of  the  form 

9.<h(T)>  vanish  by  virtue  of  Eq.  (4-20).   Therefore,  in 
J 
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particular, 


<3  . 

J 


1+T 


>  =  0  . 


C6--U ) 


From  this  follow  the  identities 

2  -i 


2(3.t) 

J 

(1+T)JJ 


>  =  < 


3?T 


(1  +  T) 


(6-5) 


and 


3,T 

J 


L(1+t) 


>  =  0 


(6-6) 


Thus,  with  Eqs.  (6-5)  and  (6-6),    the  average  of  Eq: 
and  (6-3)  are 
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1    ) 


(6-8) 
cont . 
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R' 


k     <IJL_i>   =   - 


1+T 


k  2R"       R'    . 


(6-8) 


R' 


R  R' 


If   Eqs.    (5-19)    which   give   expressions    for   T° 0    and   T   fl 

Pd     Eos      (6-7)    and    (6-8)    become    (with  the   aid   of 
are  used ,    tqs .    v."    i  * 

Eq.    (5-22))  ^    _^  [fs.il    1 

k<D(1   +   ctu2)> 


+  3RJ_  < 
+      R 


'(x+x)2 


X+t 


5      -G 


<v> 
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r2sin26      L(1+t)    J 


3k        3RJ. 


»-*«i*i' 
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3R1     , 
J£  <Du2(2a-l)    -    3aD>    +   — 
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(1+T 


>     - 


IT 


(l+x) 


(6-9) 

"1 


+   < 


R2      l 


k_        2R"       R_l 
R2  R         R2 


(6-10) 


These 


two  equations  are 


the  desired  dynamic  equations. 
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B.   CONSTRUCTION  OF  THE  SPECIAL  MOMENT  EQUATIONS 

As  stated  in  Ch.  II  ,  it  is  necessary  to  determine  the 
functional  form  of  the  extra  source  terms  of  Eqs.  (6-9) 
and  (6-10)  in  order  to  solve  for  the  average  quantities  R, 


0  >,  and  <TJ .>. 

The  basic  idea  is  to  construct  additional  equations 


for  the  unknown  correlations.   This  is  done  by  first  multi- 
plying the  unaveraged  Eqs.  (6-2)  and  (6-3)  or  combinations 
of  these  by  some  "special  functions"  which  involve  powers 
or  derivatives  of  t  and  then  averaging  them.   As  might  be 
expected,  the  nonlinearities  of  the  unaveraged  equations 
cause  more  unknown  correlations  to  be  formed  which  requires 
the  formation  of  more  equations.   Although  this  process  can 
be  carried  out  ad  infinitum,  a  special  truncation  of  this 
process  is  required  which  imposes  an  additional  restriction 
upon  the  model.   The  truncation  that  will  be  used  is  the 
higher-order-moment  discard  approximation  that  is  discussed 
in  Ch.  Ill  in  connection  with  the  fluid  turbulence  hierarchy 
This  means  that  the  higher  order  moment  that  results  from 
the  nonlinearity  is  discarded  because  of  smallness.   Hence, 
only  the  first  "generation"  of  the  moment  hierarchy  is 
needed  and  it  turns  out  that  three  special  equations  are 
sufficient  (see  Ch.  VH). 

The  two  special  combinations  of  the  unaveraged  dynamic 
Eqs.  (6-2)  and  (6-3)  which  are  used  as  the  basis  for  develop- 
ing these  three  equations  are 
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k(T°0  +  TJ  .)  =  G°0  +  G« 
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(6-12) 


If  Eq.  (6-11)  is  multiplied  by  the  "special  function" 
l+T  and  averaged,  the  first  moment  equation  is  formed.   The 
resulting  equation  can  be  simplified  by  again  using  Eq . 
(4-20)  which  implies  that 

82<h(T)>  =  0  . 
J 


If  h(i)  is  taken  to  be  ln(l+x),  then  the  following  identity 
results 


± 

H+tJ 


>  =  < 


(8,t) 
(1  +  t) 


(6-13) 


Thus,  the  simplified  first  moment  equation  is 
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If  Eq.  (6-12)  is  multiplied  by  T  and  then  8oT,  the 
second  and  third  moment  equations  are  produced. 
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Again,  using  equation  set  (5-19),  these  three  equations 
become 
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VII.   TRUNCATION  OF  THE  MOMENT  EQUATION  HIERARCHY 

In  both  Eqs.  (6-9)  and  (6-10)  it  can  be  seen  that  the 
source  terms  contain  (1+T)n  in  their  denominators.   Because 
(1+x)   can  be  expanded  in  a  power  series  in  x,  the  source 
terms  can  be  expressed  as  an  infinite  sum  of  moments  of 
various  orders  in  t  and  its  derivatives.   As  mentioned 
before,  an  equation  for  each  of  these  moments  can  be  gener- 
ated by  multiplying  the  unaveraged  Eqs.  (6-2)  and  (6-3) 
by  an  appropriate  quantity  and  averaging.   The  nonlinearities 
of  these  equations  create  a  dependence  of  the  desired  moment 
on  higher  order  moments.   Because  this  occurs  for  all  moments, 
a  method  is  needed  to  truncate  the  system  in  such  a  way 
that  only  a  finite  number  of  moments  and  required  equations 
in  the  hierarchy  are  left.   In  Ch.  Ill  three  possible  trun- 
cation schemes  are  presented.   In  this  work,  the  higher- 
order-moment  discard  method  is  chosen  by  assuming  that 
most  of  the  probability  mass  of  the  random  variable  x(x  ,£) 
is  concentrated  around  x  =  0,  so  that  the  higher  order 
moments  can  be  neglected.   This  obviously  limits  the  appli- 
cability of  the  theory  to  models  of  small  nonunif ormity 
for  which  |t|<<1  (see  Ch.  IV-B) . 

However,  care  must  be  exercised  in  affecting  this  trun- 
cation.  It  is  not  correct  to  just  replace  (1+x)   by  one  in 
each  of  the  denominators  of  the  above  mentioned  equations 
just  because  this  is  the  first  term  of  the  expansion. 
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can  be  seen  by  an  example.   Consider  the  typical  term 
This  term  can  be  expanded  as  follows 


>  =  Oqt>  -  <(30t)(t)>  +  .  .  . 


fa  T1 
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(7-1) 


Since  O  x>  =  9  <t>  =  0, 
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(7-2) 


However,  if  (1+t)  were  falsely  replaced  by  one,  the  follow- 
ing incorrect  conclusion  would  result  which  disagrees  with 
the  correct  Eq .  (7-2) 
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O  T> 


(7-3) 


Hence,  an  expansion  must  be  made  for  every  term  and  that 
lowest  order  moment  term  which  does  not  vanish  by  virtue 
of  having- a  zero  average  is  assumed  to  be  a  "good"  approxi- 
mation to  the  term  itself. 

With  this  "weak  fluctuation"  assumption,  Eqs.  (6-9) 
(6-10),  (6-17),  (6-18),  and  (6-19)  become 
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In  the  limit  of  zero  fluctuations,  Eqs .  (7-4)  and  (7-5) 
become  the  usual  two  cosmological  equations  for  p  and  R 
that  are  given  by  Adler  [1].   Thus,  it  is  natural  to  look 
upon   these  equations  as  dynamic  equations  for  p  and  R. 
The  source  terms  which  appear  in  these  equations  contain 

2  2 

the  three  functions  <t2>  ,  <(3„t)  >,  and  {<(8  t)  > 


+  ?T  <(V)2>  +  rasing  <(V}  >}'   ThUS'  0nly  ^^ 
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equations  are  needed  to  determine  these  variables  which 
justifies  the  statement  made  in  Part  B  of  the  preceeding 
chapter  that  only  three  equations  of  the  hierarchy  need 
to  be  calculated.   These  three  equations  are  Eqs .  (7-6), 
(7-7),  and  (7-8).   Although  the  above  functions  <x2>, 

2  2       -,  2  -,  2 

<(8,t)  >,  and  {<(9  t  )  >  +  ±y  <(3ax)  >  +   ,  .  ,a    <Oax)  >} 
os      L    r       r2     6       r2sm28    <J>     ' 

appear  in  these  equations,  additional  moments  also  occur 
which  relate  the  density  D  to  x  and  D  to  9Qx.   It  is  shown 
in  the  succeeding  chapters  that  these  moments  can  be  related 
to  the  existing  moment  <x2>.   Thus,  Eqs.  (7-k)    to  (7-8) 
are  soluble  among  themselves.   Because  Eqs.  (7-6)  to  (7-8) 
govern  the  dependence  of  the  source  moments  on  R,  they 
are  called  growth  equations  to  convey  the  idea  that  they 
determine  the  development  or  growth  of  a  perturbation  x  with 
increasing  R. 
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VIII.   REVIEW  OF  THE  RESTRICTIVE  ASSUMPTIONS 

Before  proceeding  further,  it  might  be  helpful  to  reviev; 
the  restrictive  assumptions  that  are  used  to  obtain  Eqs. 
(7-^)  through  (7-8).   The  first  is  the  very  restrictive 
form  of  the  metric  tensor.   This  form  excludes  many  types  of 
gravitational  wave  phenomena.   All  that  is  seen  by  an  observ- 
er would  be  a  local  "expansion"  or  "contraction"  of  space 
because,  as  stated  in  Ch.  IV,  it  would  always  be  possible 
to  construct  an  expanding  or  contracting  Robertson-Walker 
geometry  that  is  "tangent"  to  the  given  space  at  any  given 
point  in  space-time.   The  exclusion  of  such  waves  must 
certainly  exclude  any  hope  of  using  this  model  to  describe, 
in  a  realistic  way,  the  violent,  multiply  connected  space 
that  Wheeler  envisions  must  typify  the  quantum  fluctuations 
of  space-time. 

The  second  and  third  restrictions  are  the  restrictions 
concerning  the  Robertson-Walker  form  of  the  average  metric 
and  the  spatial  stationarity  of  the  fluctuating  part.   These 
two  assumptions  rely  more  on  intuition  than  upon  scientific 
rigor.   They  imply  that  there  exists  some  sort  of  spatial 
uniformity  in  the  statistics  of  the  assumed  fluctuating 
geometry.   However,  it  is  to  be  understood  that  a  "new" 
cosmological  principle  is  not  being  seriously  proposed  here 
by  the  use  of  these  two  assumptions.   That  these  assumptions 
are  used  at  all  reflects  a  need  to  simplify  the  mathemat 
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of  the  moment  equations  in  some  reasonable  way  by  appealing 
to  some  plausible  symmetry  arguments.   Thus,  these  two  re- 
strictions should  be  viewed  only  as  steps  in  this  particular 
problem  whose  purpose  is  to  demonstrate  a  technique. 

Finally,  the  fourth  restriction   is  the  assumption 
that  the  higher  order  moments  are  negligible  owing  to  an 
assumed  smallness  in  the  anisotropy  and  inhomogeneity .   This 
assumption  is  a  restriction  upon  the  applicability  of  the 
model.   That  is,  in  those  solution  domains  for  which  the 
magnitude  of  t  is  comparable  or  larger  than  one,  the  use  of 
the  model  is  unwarrantable .   If  large  amplitude  fluctuations 
are  to  be  considered,  then  perhaps  a  truncation  based  upon 
a  phenomenological  or  a  cumulant-discard  approximation  would 
be  useful.   However,  in  the  case  of  the  latter,  it  might 
be  assumed  that  the  probability  distribution  of  the  random 
variable  x  is  normal.   Then,  the  higher  order  moments  can 
be  shown  to  be  products  of  the  lower  order  ones  [18].   Such 
an  assumption  is  made  in  fluid  turbulence  theory  [3l«   But, 
in  the  case  of  this  model,  x   must  be  greater  than  or  equal 
to  minus  one  as  given  by  Eq.  (4-15).   This  creates  difficulty 
with  a  normal  distribution  where  the  probability  of  x  being 
+T  is  equal  to  the  probability  of  it  being  -T  where  T  is  an 
arbitrary  positive  number. 
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IX.   SPECIAL  UNIVERSE  MODELS 

Prom  this  point  on  it  is  convenient  to  consider  two 
types  of  models  which  correspond  to  two  values  of  a.      The 
first,  for  which  a  is  equal  to  one  third,  corresponds  to  a 
universe  made  out  of  electromagnetic  radiation  while  the 

second  corresponds  to  a  universe  which  consists  of  cold 

p 

matter  for  which  a  is  zero.   (a  is  equal  to  —tit  •   See  Eq. 

^       c  2D         ^ 

(5-22).) 

It  is  not  practical  to  treat  these  models  simultaneously 
by  leaving  a  arbitrary.   The  cold  dust  model  can  be  solved 
exactly  owing  to  the  fact  that  all  matter  moves  freely  along 
geodesies  without  being  accelerated  by  any  pressure  gradi- 
ents (pressure  is  zero  everywhere).   But  the  radiation  model 
does  involve  pressure  gradients  which  do  not  allow  the 
elements  of  the  radiation  fluid  to  move  freely  in  space   [9] 
Therefore,  these  two  types  of  models  are  sufficiently  differ- 
ent from  one  another  to  merit  separate  consideration. 


A.   RADIATION  FILLED  UNIVERSE  MODEL 

For  the  radiation  filled  universe,  Eqs .  (7-^),  (7-5), 
(7-6),  (7-7),  and  (7-8)  become 
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where 


and 


p  e  <D(1  +  |-)> 


d  =   D(l  +  H_)  _.  p 


(9-6) 


(9-7) 


It  should  be  noted  that  P  is  equal  to  the  average  of  T  Q  . 
Adler  [1]  points  out  that  T   can  be  interpreted  as  the 
total  mass  density  with  respect  to  the  given  coordinate 
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system  and  therefore  includes  relativistic  corrections  to 
the  proper  mass  density. 

1.   Elimination  of  the  Spatial  Derivatives 

Nov;  consider  the  terms  which  occur  in  Eqs.  (9-1) 

2 

and  (9-2)  of  the  form  <(9  t)  >.   It  should  be  noted  that 
only  9Q  and  9   are  partial  derivatives  with  respect  to  a 
length  (in  the  time  and  r  directions).   The  remaining  par- 
tial derivatives  9Q  and  9  .  are  partial  derivatives  with 
respect  to  the  coordinates  9  and  <f>  which  are  not  lengths 
in  themselves.   An  inspection  of  the  line  element  of  Eq. 
(4-2)  suggests  that  9Q,  e~    9  ,  e~    1/r  9Q,  and  e~ 
l/rsin8  9,,  all  of  which  have  dimensions  of  inverse  lengths, 
play  the  same  roles  in  the  x°,r,6,  and  <J>  directions  respec- 
tively.  Because  disturbances  of  the  metric  propagate  with 
the  speed  of  light  [11,  32],  it  would  be  expected  that  the 
characteristic  inverse  lengths  associated  with  these  quanti- 
ties are  comparable  in  magnitude.   Thus,  it  seems  reasonable 
to  expect  that  the  spatial  terms 

e"G<0  t)2>,  i-e-G<(3  x)2>(  — S _2<(3  t)  \ 

r  (r  sin9) 

could  be  expressed  as  mere  repetitions  of  the  temporal  term 

2 

<(9  x)  >  multiplied  by  a  coefficient  of  order  one. 

Such  a  fact  would  simplify  Eqs.  (9-D    and  (9-2) 
because  the  terms  involving  the  spatial  pertial  derivatives 
would  be  changed  into  the  already  existing  terms  involving 
the  time  derivative. 
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This  consolidation  of  the.  terms  is  performed  in 
Appendix  F.   Equation  (9-3)  is  used  to  eliminate  the  spatial 
partial  derivative  terms  from  Eqs .  (9-1)  and  (9-2).   An 
appropriate  renormalization  is  then  performed  to  bring  the 
right  hand  sides  of  the  resulting  equations  back  to  "canoni- 
cal form"  as  they  appear  in  the  right  hand  sides  of  Eqs. 
(9-1)  and  (9-2).   Equation  (9-*0  is  used  to  change  the  form 
of  the  "renormalized"  equations  into  a  more  convenient  form. 
This  last  step  is  not  part  of  the  consolidation  itself 
since  it  just  allows  the  resulting  equations  to  be  presented 
in  a  more  desirable  form.   The  transformed  pair  of  Eqs. 
(9-D  and  (9-2)  which  are  void  of  spatial  partial  deriva- 
tives are 

kP  +  f  K<id>  +  }   <(30t)2>  +  |  j^-  90<x2> 
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2.   Inclusion  of  the  Dynamics  for  the  Radiation  Field 
Because  Eq.  (9-3)  has  been  used  to  remove  the  un- 
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from  Eqs.  (9-1)  and  (9-2),  there  exists  now  only  four 
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equations  In  the  system  namely,  Eqs .  (9-4),  (9-50,  (9-8), 
and  (9-9).   But  this  system  of  equations  is  insoluble  because 
they  represent  four  coupled  nonlinear  differential  equations 
in  six  unknown  functions  of  the  time  coordinate  x°.   These 

2 

unknown  functions  are  p,  Rs  <t2>,   <0  t)  >,  <xd>,  and  <0  x)d>. 

However,  there  is  a  way  to  express  <xd>  and"  <0 -i)d> 
in  terms  of  the  other  four  unknown  functions.   It  is  not 
unreasonable  to  expect  that  there  should  exist  some  rela- 
tionship between  the  density  fluctuation  d(x  )  and  geometri- 
cal fluctuation  t(x^)  because  in  some  sense  the  geometry 
fluctuations  "drive"  the  density  fluctuations.   But  it  must 
be  remembered  that  the  radiation  fluid  possesses  internal 
pressure  forces  which  push  the  fluid  elements  about  as 
well  as  do  work  against  any  expansions  or  contractions  of 
space  which  might  occur  with  amplitude  t  [9].   Thus,  the 
problem  of  determining  this  relationship  between  d  and  t  is 
in  general  very  complicated.   The  correct  way  to  do  this 
is  to  use  the  equations  of  motion  for  the  radiation  fluid 
which  are  given  in  tensor  form  by  [1] 

T^v;  v  =  0  (9-10) 

where 

T^v:  v  =  8  T^v  +  rV   Tya  +  ry  Tav  . 
'      v       va        av 

Indeed,  the  relationship  between  d  and  t  is  approximated 
in  Appendix  I  via  these  equations  and  is  expressed  in  terms 
of  the  correlations  <id>  and  <(30i)d>. 
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An  instructive  yet  simpler  way  to  determine  the  rela- 
tionship between  d  and  t  is  the  following.   Consider  a  small 
3-volume  element  V  at  rest  with  respect  to  the  coordinates 
(r,  6,  <J>).   Let  M  be  the  total  relativistic  mass  content 
of  this  3-volume.   The  total  relativistic  mass  density  D 
is  thus  given  by 

D  =  |    .  (9-1D 

A  small  change  AD  of  D  can  be  in  general  a  result  of  a  small 
change  AM  of  M  and  a  small  change  AV  of  V.   That  is, 

AD  =  AM  _  MAV   ^  (9_12) 

V     V2 

Now  the  quantity  /^g"  d^x  is  known  to  be  a  scalar 
invariant  where  g  is  the  determinent  of  the  matrix  g 
[1,  32].   In  particular,  from  special  relativity,  this  quan- 
tity is  the  product  of  the  proper  time  interval  and  the 
proper  volume.   Because  x°  has  already  been  identified  as 
the  time  coordinate  (see  Ch.  IV),  the  volume  V  must  be 

V  =  JZg   d3x  .  (9-13) 
In  terms  of  the  coordinates  r,  6,  cf> , 

V  =  /Zg   drdOdcj)  .  (9-1*0 
The  use  of  Eq.  (4-8)  gives 

V  =  e3M/2(dr)(rde)(r  sined^).  (9-15) 
But     eM  =  eG(l+x)  (9-16) 
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so  that  V  =  V(1+t)3/2  (9-17) 

where   V  =  e3G/2 (dr ) (rde ) (r  sin0d<j))  . 

The  change  of  mass  AM  can  occur  because  of  the  work 
done  by  the  pressure  against  the  change  of  volume  and  by 
the  flux  of  mass  out  of  the  volume.   If  this  flux  is  tempo- 
rarily ignored  (i.e.  all  mass  moves  along  geodesies),  then 
AM  can  be  expressed  as 

AM  =  -  —  AV  (9-18) 

c2 

where  P  is  the  total  relativistic  pressure.   The  relation 

P  =  |  c2D  (9-19) 

implies  that 

AM  =  -  i  M  ^  .  (9-20) 

Then  Eq.  (9-12)  becomes,  using  Eq .  (9-20) 

AD  =  -  |  D  f-    .  (9-21) 

But  Eq.  (9-17)  with   small  becomes 

V  =  V  +  |  V x  .  (9-22) 

If  AV  is  replaced  by  |  Vt ,  V  by  V,  AD  by  d,  and  D  by  p  in 
Eq.  (9-21),  then  this  equation  becomes 

d  =  -  2Pt    .  '  (9-23) 

Hence,    the   correlation   function   <xd>    and    <O0i)d>   become 
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<id>  =  -  2p<x2> 
and 

<O0x)d>  =  ~  P80<T2>   .  (9-24) 

In  Appendix  I,  the  effect  of  "mass  leakage"  from  the  volume 
is  approximated  by  means  of  two  variable  parameters  X  and  6. 
It  is  shown  that  a  reasonable  form  for  the  correlations 
<xd>  and  <(8  x)d>  are 

<xd>  =  -  2Ap<i2> 

<O0T)d>  =  -  ep30<x2>  .  (9-25) 

The  coefficients  A  and  $  may  be  thought  of  as  representing 
a  zeroth  approximation  of  those  functions  of  x°,  which 
correct  the  zero  leakage  correlation  Eqs .  (9-24).   These 
coefficients,  then,  become  open  parameters  which  character- 
ize a  given  model.   Certainly,  if  mass  leakage  were  not 
important,  then  A  and  3  would  both  be  one. 

The  substitution  of  Eqs.  (9-25)  into  the  dynamic 
Eqs.  (9-8)  and  (9-9)  and  into  the  growth  Eqs.  (9-4)  and 
(9-5)  gives 

kP  -  |  AkP<t2>  +  }   <O0t)  >  +  |  jp  80<t2> 
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8o<x2>  +^TVt2>  +|akP<t2>  =  2<(8ot)2>    (9-28) 

2    Mr '      2      M 
30<Oot)  >  +  —  <(30t)  >  "  "  J  6kP80<t2>.     (9-29) 

3-   Solution  of  the  Moment  Equations  which  Govern  Growth 
of  the  Fluctuations 
The  four  Eqs.  (9-26)  through  (9-29)  now  contain  the 

2 

four  variables  p,  R,  <t2>,  and  <(8qt)  >  and  the  system  is 
now  closed.   However,  the  solution  of  these  four  equations 
is  not  easy  because  they  are  nonlinear  and  nonlinear  prob- 
lems do  not  in  general  yield  easily  to  analysis.   But  it 
is  possible  to  take  advantage  of  the  requirement  that  the 
moments  which  appear  in  Eqs.  (9-26)  and  (9-27)  be   small 
compared  to  xp.   That  is,  the  additional  stress  terms  cause 
a  small  perturbation  of  the  usual,  fluctuation  free  cosmol- 
ogical  solutions  for  p  and  R  which  are  well  known.   These 
solutions,  derived  in  Appendix  E,  are 

3a 
kP  = £  (9-30) 


R1* 


and 


R.2  =  i  _  k  (9-3D 

R2 

where  a  is  a  constant  of  integration  and  k  is  the  curvature 
r 

characteristic  and  may  take  one  of  three  values  -1,  0,  +1 
(see  Ch.  IV). 

Nonlinear  problems  can  in  many  cases  be  simplified 
by  the  replacement  of  a  variable  with  a  very  close  approxi- 
mation to  that  variable  [30].   If  such  a  replacement  is 
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made  for  R'  and  p  in  Eqs .  (9-28)  and  (9-29)  with  the  fluc- 
tuation free  solutions  Eqs.  (9-30)  and  (9-3D  and  if  all 
derivatives  with  respect  to  x°  are  transformed  into  deriva- 
tives with  respect  to  R  by  the  chain  rule  for  differentia- 
tion, these  two  equations  become  soluble  between  themselves. 
However,  as  can  be  seen  from  Eq .  (9-31),  the  effect 

2 

of  k  upon  R1   is  not  significant  unless  R  is  large  enough 
so  that  ar/R2  is  of  the  same  order  of  magnitude  as  k.   It 
is  well  understood  that  that  value  of  R  for  which  the  effect 
of  k  is  important  to  the  dynamics  is  the  current  value  of 
R.   For  this  reason,  much  research  in  observational  cosmol- 
ogy Is  devoted  to  the  determination  of  k  on  the  basis  of 
data  on  current  expansion  characteristics  [25].   But  the 
applicability  of  the  radiation  model  is  only  valid  up  to 
the  time  when  the  radiation  density  is  equal  to  the  matter 
density  for  which  the  associated  R-value  is  much  less  than 
the  current  value  of  R  by  several  orders  of  magnitude  [8]. 
Thus,  if  the  R-domain  of  the  solutions  <t2>(R)  and 
<(90t)2>(R)  is  restricted  to  values  of  R  sufficiently  less 
than  the  current  value  of  R  so  that  ar/R2»  |k| ,  then  the 

9 

following  substitutions  for  R'   and  p  in  Eqs.  (9-28)  and 
(9-29)  are  valid. 

3a 
Kp  «  -  -£  (9-32) 

R" 


and 


R.2  =  !L  (9-33) 

R2 
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It  should  be  noted  that  the  neglect  of  k  in  cosmol- 
oglcal  problems  is  not  new.   k  is  neglected,  for  example, 
by  Hoyle  and  others  in  their  expression  for  R'   for  exactly 
the  same  reason  that  it  is  neglected  here  [10]. 

The  growth  equations  become,  using  the  above  approx- 
imations 

R28^<T2>  +  R3r<t2>  -  8x<x2>  =  |^-  <O0x)2>     (9-3*0 


R 


<(30t)  > 


+  i 


<O0t)  > 


=  ^a.<x2>. 


R 


<♦  R 


(9-35) 


These  equations  are  solved  in  Appendix  J  for  <x  > 


and  <(3  t)  >.   The  solutions  are 


<x2>  =  F 


+  V3  Ri^  +  ^Z3r^ 


1/3    ^3 


H 


(9-36) 


—1 .  -  4XF1/3R   +  I^R1  +  D1/3R^ 


(9-37) 


where   i  =  -  4  +  2/20  +  3) 
j  =  _  l\   _  2/2(X  +  3) 

and  I. /q  ,  D,  ,,  ,  Fi/-z   are  arbitrary  constants  of  integration, 

4  .   Special  Submodels 

a.   Discussion 

In  general,  a  given  combination  of  the  parameters 

T      D      and  F  ,      specify  a  nonuniform  model  which 
1l/3  »   1/3  '      1/3  y 

is  difficult  to  treat  analytically.   Yet,  there  are  special 
values  of  these  parameters  which  designate  models 
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relatively  easy  to  treat  and  which  are  interesting  .  These 
are  the  three  models  for  which  only  one  of  the  three  arbi- 
trary parameters  are  non  zero.   That  is, 

<(90t)    >  h 

(1)  <t2>   =   F1/3    ;    — — =    -    4AF1/3R-4  (9-38) 

2 
In   /q  -t  j.Ji  <(80T)      >  , 

(2)  <.t'».^R1+«    ;-4-      --I^1  (9-39) 


2 

(3)    .^fR^'    ;   — ^_-=D1/3rJ 


(9-40) 


Actually,  not  much  generality  is  lost  by  considering  the 
above  three  models  separately.   The  I  ,   ,  D  ,   ,  and  F  , 
terms  of  the  general  Eqs.  (9-36)  and  (9-37)  are  functions 
of  R  to  a  power.   In  general,  however,  these  powers  are 
quite  different  from  one  another  (i.e.  try  X,  3  equal  to 
one  which  are  the  zero  leakage  values  for  A  and  3).   Because 
R  changes  by  many  orders  of  magnitudes  during  expansions 
and  contractions,  only  one  of  the  terms  would  dominate 
anyway  (except  for  those  "particular"  values  of  R  for  which 
two  or  even  three  of  the  terms  would  be  of  the  same  order). 
Hence,  most  of  the  models  with  a  given  combination  of  Iwo  • 
D,  ._  ,  and  F  ,~  parameters  would  degenerate  into  one  of  the 
above  three  models  for  almost  all  R-values. 

However,  there  is  a  difficulty  with  the  F-wo" 
model.   Equation  (9-25)  says 

<xd>  =  -  2Ap<x2>  . 
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X  must  therefore  be  a  parameter  which  measures  the  response 
of  the  radiation  fluid  to  a  given  geometric  perturbation. 
It  seems  unreasonable  to  suppose  that  the  fluctuation  d  of 
the  total  mass  density  is  positive  for  positive  T.   That  is, 
it  seems  that  only  models  that  are  physically  reasonable 
are  the  ones  for  which  A  is  positive.   But,  from  Eq .  (9-38), 
only  the  negative  A  values  permit  the  obviously  positive 

2 

functions  <(8qt)  >  and  <t2>  to  be  positive  (a   is  positive 
by  virtue  of  Eq .  (9-32)).   Hence,  the  F     ,   -model  does  not 
seem  to  have  much  physical  meaning  and  is  therefore  not 
considered  in  this  work.   (If  some  positive  A  values  would 
be  acceptable  then  there  would  be  grounds  for  exploring 
the  content  of  the  P,  ,,,-model.  ) 

The  two  remaining  models,  namely  the  I,  ,~  and 
D,  /-.-models ,  are  very  different  from  one  another.   The 
I, ,_-model  represents  a  model  which  becomes  more  nonuniform 
with  increasing  R  while  the  D,  ,„-model  depicts  a  model  for 
which  the  nonunif ormity  decreases  with  increasing  R. 
Because  there  is  no  intrinsic  reason  why  positive  A  values 
are  not  acceptable,  both  of  these  models  can  be  considered. 
(I,,...  is  used  to  denote  the  arbitrary  constant  of  an 
increasing  type  perturbation  field  of  the  metric  for  a 
radiation,  a  equal  1/3  universe  and  D-wo  is  used  to  denote 
the  decreasing  type.   The  same  nomenclature  will  be  used 
to  specify  the  increasing  and  decreasing  type  perturbations 
of  the  dust,  a  equal  zero  universe  but  with  a  zero  subscript.) 
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It  is  of  interest  to  see  how  the   fluctuation 
stresses  which  correspond  to  both  the  I,  .~   and  D  .   -type 
perturbations  act  upon  the  average  total  mass  density  p 
and  the  average  scale  length  R.   This  can  be  seen  if  the 
stress  terms  which  appear  with  the  average  density  p  in 
Eqs.  (9-26)  and  (9-27)  are  expressed  in  a  form  which  cor- 
responds to  the  functional  form  of  the  I,  ,^   and  D,  ,_-type 
fluctuations  given  by  Eqs.  (9-39)  and  (9-40).   Equations 
(9-26)  and  (9-27)  become  (with  |y  <j2>    neglected  to 
kp<t2>  for  reasons  already  stated) 

(I, /o-f luctuations) 

™  4.  r(i         W1    [60X  +  243  +  6/2(A+3)] 
Kp  +  L{1^,^)H      j2 

3k   3R' 


R2     R2 


KP    nfT    ^pi  [12X  +  83  -  2/2U+3)] 
3~  ~   ai/3;  72 


2 

k     2R"   R1 


Rz     R    R' 


(D,  /r,-fluctuations) 


n(T>        ^pJ    [60A    +   243    -   6/2(A+3)] 
Kp    +   C(D,  ^-JR"    j2 


2 


R2  R2 


K 


-  §■  p    -   C(D,  ,.)R 


j    [12A    +    83    +    2/2(A+3)] 


(9-^D 


(9-^2) 


3i£  _   3R1_  (9-43) 


2 

k_  _     2R_^'  _   R±_  (9-44) 

R2  R  R2 
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where   0(1^)  =  -- ^-^  and  0(0^)  =   -^ 


i  =  -  4  +  2/2(X+3)   and  j  =  -  4  -   2/2(X+3)  ,   (9-45) 

It  should  be  pointed  out  that  of  the  three 
parameters  C(A),  X,  and  3  of  Eq.  (9-45)  (v;here  A  represents 
^1/1   °r  Dl/^  '  only  C  lends  itself  to  physical  measurement. 
It  must  not  be  forgotten  that  C.  is  proportional  to  a  para- 
meter v/hich  is  the  result  of  an  integration  and  therefore 
is  an  arbitrary  constant  which  must  be  fixed  by  some  obser- 
vation.  But  X   and  3  are  not  the  result  of  any  integration. 
They  are  variable  parameters  which,  in  a  sense,  measure  the 
mathematical  uncertainty  between  two  pairs  of  correlation 
functions  (see  Eqs.  (9-25)).   If  it  were  known  how  to 
exactly  determine  these  relations,  then  3  and  X  would  be 
known  numbers  and  hence  be  of  nonvariable  status.   As  it 
is,  they  have  to  exist  as  variable  parameters.   (In  the 
dust  model  considered  later,  the  dynamics  of  the  matter 
are    determined  exactly.   Hence,  the  correlation  functions 
can  be  found  precisely  and  X  and  3  do  in  fact,  for  that 
case,  become  known  numbers.) 

The  solution  of  the  dynamic  Eqs.  (9-^D  to 
(9-44)  are  given  in  the  next  section. 
b.   Classification 

The  question  naturally  arises:  What  is  the 
effect  of  a  given  combination  of  fluctuation  parameters 
C(A),  X,  and  3  upon  p  and  R?   It  has  already  been  pointed 
out  in  Section  3  of  this  chapter  that  the  effect  of  the 
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fluctuation  stresses  upon  the  average  density  p  and  scale 
length  R  is  small.   Thus,  the  actual  solutions  are  small 
perturbations  of  the  fluctuation  free  solutions  and  it 
tnerefore  seems  appropriate  to  compare  the  solutions  of 
Eqs.  (9-^D  through  (9-W  for  p  and  R  with  these  known 
fluctuation  free  solutions. 

This  can  be  done  in  the  following  manner:  There 
are  only  four  distinct  ways  that  the  fluctuation  stresses 
can  perturb  the  solutions  for  p  and  R  from  their  fluctuation 
free  solutions.   That  is,  at  a  given  time  x°,  the  density 
p  and  the  scale  length  R  can  each  be  perturbed  above  or 
belov;  their  fluctuation  free  values.   Let  each  of  these 
four  type  of  perturbations  be  designated  by  the  letters  K, 
L,  M,  and  N  in  the  manner  shown  in  Fig.  1  for  the  increasing 
I,  ,_-f luctuations  and  Fig.  2  for  the  decreasing  D,  /-.-fluctu- 
ations.  In  both  cases,  a  given  letter  indicates  whether 
the  p  and  R  are  greater  or  smaller  than  the  fluctuation 
free  solutions  for  p  and  R  respectively.   For  example,  a 
model  for  which  p(x°)  and  R(x°)  are  both  greater  than  their 
respective  unperturbed  values  p*(x°)  and  R*  would  be  clas- 
sified as  type  M. 

The  manner  in  which  the  solutions  p  and  R  are 
perturbed  from  their  fluctuation  free  values  which  are 
specified  by  the  above  simple  classification  scheme  is 
determined  solely  by  the  two  parameters  X  and  3.   The  mag- 
nitude by  which  a  perturbation  of  a  given  type  occurs  is 
determined  by  the  parameter  C(A).   Thus,  X  and  3  determine 
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the  qualitative  behavior  of  the  model  while  C(A)  determines 
its  quantitative  behavior. 

In  order  to  determine  the  classification  of  a 
model  which  is  characterized  by  the  three  parameters  C(A), 
X,  and  3  where  A  refers  to  either  the  increasing  I,  ,,-type 
perturbation  or  the  decreasing  D,/--type,  the  solutions  for 

p  and  R  are  needed.   In  Appendix  L  the  solutions  for  <p 

2 
and  R'   are  given  in  terms  of  some  special  parameters  f . , 

bA,  and  m.  which  are  functions  of  only  the  three  parameters 

C(A),  A,  and  6.   These  solutions  are 

~H  mA 

Kp  -  -  3arR   -  f AR  A  (9-^6) 

2       -2       mA+2 
R'   =  arR   +  bAR  A     .  (9-217) 

The  unperturbed  solutions  p*  and  R*  are  those 
for  which  f.  and  b.  are  both  equal  to  zero.   It  is  easily 
seen  that  the  direction  of  displacement  of  p  from  p*  is 
determined  solely  by  the  sign  of  f..   The  sign  of  the  para- 

2 

meter  b.,  on  the  other  hand,  tells  whether  R'   is  larger 

2K 

or  smaller  than  its  unperturbed  value  (RT)   for  all  values 
of  R.   If,  in  the  case  of  the  I  ,^-type  model,  a  forward 

time  integration  of  Eq.  (9-^7)  for  various  values  of  bT 

o   1/3 

is  made  from  R  equal  zero,  a  family  of  solutions  R(x  ) 

can  be  generated.   The  displacement  of  any  one  of  these 

2      2  * 

solutions  from  R*  depends  on  the  ratio  R'  /$. '  )  .   For  example, 
if  this  ratio  is  greater  than  one  for  all  R-values,  then  R 
is  necessarily  perturbed  above  R*  because  all  solutions  of 
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the  family  are  "pinned"  to  the  same  initial  value.   Because 

the  size  of  this  ratio  compared  to  one  is  determined  by  bT 

Xl/3' 

it  follows  that  bT    also  determines  the  direction  of  the 

1l/3 

displacement  of  R  from  R* .   It  can  be  shown  that  the  sign 
of  t>D-j/-3  also  determines  the  displacement  of  the  decreasing 
D,  ,-,  solution  of  R  from  R*  if  the  integration  is  begun 
at  a  sufficiently  large  R  and  performed  backwards  in  time. 
In  Appendix  L,  the  signs  of  f.  and  b.  (and  hence  the  classi- 
fication letters)  are  determined  for  all  values  of  X  and  3 
for  both  the  D,  ,_  and  I  ,„  type  fluctuations.   The  results 
are  graphically  displayed  in  Figs.  3  and  4. 

It  might  be  asked  how  a  given  combination  of   X 
and  $  affect  the  exponent  of  R  in  the  fluctuation  stress 
term.   The  dependence  of  this  exponent  is  presented  in  Fig. 
5  as  a  plot  of  the  lines  of  constant  j  and  i  in  X-8  space. 

This  concludes  the  basic  discussion  and  calcu- 
lations for  the  radiation  model.   In  the  next  section,  a 
similar  treatment  is  made  for  the  dust  model.   The  discus- 
sion and  comparison  of  these  two  models  are  given  in 
Chapter  X. 

B.   DUST  FILLED  UNIVERSE  MODEL 

Some  of  the  calculations  and  discussions  in  this  section 
are  similar  to  ones  of  the  previous  section  which  deals 
with  the  radiation  filled  universe  for  which  a  is  equal  to 
one  third.   Here  the  model  for  which  a  is  equal  to  zero  is 
treated.   It  is  assumed  that  the  reader  is  already  familiar 
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with  the  previous  section  so  that  otherwise  redundant  expla- 
nations may  be  omitted. 

For  the  pressure  free  universe,  Eqs.  (7-4),  (7-5),  (7-6), 
(7-7),  and  (7-8)  become 


3  R' 

kp   -  j  r- 


a0<i2>  +  f  <(80t) 


5    0-G, 
If  e 


<(3t,t)    >    +   —  <(3dt)2>    + —  <(3Ax)    > 

r  r2  e  r2sin26  * 


3k<T2>   =   _   3k  _   3RJ. 
R2  R2  R2 


(9-^8) 


||i  3o<t2>   _   |  32<t2>    +   3   <(3oT)2> 


5      -G, 

la  e 


i 


<o  T)2>  +  \  <o9t)2>  +  -t-7T"<(V)  * 

1  r2  D  r2sin2fi 


k       .    2.  k  2RM        R' 


+     <TZ>      =     - 

R2  R2  R  R2 


(M9) 


5      -G, 
?  e 


<Ot,t)    >    +   —  <(3at)    >    + <(3at)    > 

r  r2  6  r2sin20  * 


5              5     ^    -          15k        15R'          15R"         ,Q   c-nx 
£  Kp   -  4  K<id>    -  -^—  -  -  (9-50; 

22  R2  R2  R 


K<xd>    -   I  32<t2>    - 


2    uo 


R      uo 


3  „<t'>    = 


-    3<(3qt)    > 


9    <Onx)2>    +   %^-  <(3fti)2>    =   |  «(30T)d> 


0       x      0 


R 


(9-5D 
(9-52) 


where 


p    =    <D> 


and 


d   E    D   -   p 


(9-53) 
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It  should  be  noted  that  for  a  equal  to  zero,  the  pressure 
is  identically  zero.   Hence,  there  are  no  pressure  gradients 
available  to  accelerate  the  matter.   This  implies  that 
matter,  when  left  alone,  will  follow  geodesic  lines.   If 
these  geodesies  are  taken  to  be  geodesies  of  zero  velocity 
then  the  matter  comoves  with  the  coordinate  system  (r,6,<J>) 
and  the  above  equations  follow.   Prom  Eqs.(9-53)  it  is  seen 
that  p  is  the  average  of  the  proper  density  D  because  the 
relative  velocity  u  is  zero  in  which  case  D  is  equal  to  D 
(see  Eq.  (5-20)) . 

1.   Elimination  of  the  Spatial  Derivatives 

In  Eqs.  (9-^8)  and  (9-^9),  the  terms  which  contain 
spatial  partial  derivatives  may  be  eliminated  explicitly  in 
a  manner  analogous  to  the  way  they  are  removed  from  Eqs. 
(9-1)  and  (9-2).   This  computation  is  performed  in  Appendix 
G.   The  resulting  two  equations  are 

KP     +     |    <O0T)2>     +     f    jp     90<T2>     +     \     K<Td> 


3   k_<x2>    =   _    3k  _   3R^_  (9_5l0 

R2  R2  R2 


- 1  <^ot)2>  + 1  ir  vt2>  +  it  K<Td> 


lk,      k    2R"   R' 


2 


Tj-  -    t2>  = 


(9-55) 


R2  R2    R    R' 

2.   Inclusion  of  the  Dynamics  for  the  Matter 

In  order  to  solve  the  system  of  equations,  it  is 
necessary  to  express  <xd>  in  terms  of  existing  correlation 
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functions.   This  is  also  done  in  the  previous  section  but 

is  complicated  by  the  existence  of  nonzero  velocities. 

For  the  pressure  free  model  this  calculation  is  much  easier. 

Consider  Eq .  (9-12).   Because  the  velocities  and 
the  pressure  are  zero,  there  is  no  mechanism  to  change  the 
mass  content  of  an  expanding  volume.   Hence,  AM  must  be 
zero  and  Eq .  (9-12)  becomes 

AD--.*F.  (9-56) 

V2 

With  Eq.  (9-22)  for  small  t,  Eq .  (9-56)  becomes 
(also  with  the  usual  definitions  for  the  increments  as  in 
the  previous  section) 

d  =  -  |  px  .  (9-57) 

Hence,  <id>  and  <;8oT)d>  must  be 

<xd>  =  -  |  p<t2>  (9-58) 


and 


2 


<(30T)d>  =  \   pVt2>   •  (9-59) 


Equation  (9-57)  is  also  derived  in  Appendix  I  using 
the  conservation  equations 

TPV;  v  -  0  . 

It  should  be  noted  that  there  is  no  need  for  any 
arbitrary  parameters  A  and  3.   This  is  because  the  above 
dynamical  equations  can  be  determined  exactly. 


n 


With  Eqs.  (9-58)  and  (9-59),  the  dynamic  Eqs.  (9-5*0 
and  (9-55)  and  the  growth  Eqs.  (9-5D  and  (9-52)  become, 


Kp 


+  3-<.O.t;>  >  +  ||^  80<t2>  - 


3  H_<x2>  =  _  3k  _  3RL 
*    R2  R2     R2 


15      2 

""§"  Kp<TZ> 


-  i<(8°T)2>  +  £  jr  Vt2>  -  kkp<t2> 


(9-60) 


1  k 
R2 


<t2>  =  _  L_  _  2R_1<_  R 
R2     R     R' 


(9-61) 


9rr<t2>  + 


fR11     J  1 


R 


R 


8u<t2>  +  ^ 


'H 


R 


r  <t  > 


^<(8oT)2> 


3D<(3nT)  >  +  5"  <  (9„t)2>  = 


R 


R 


2  p9R<T2> 


(9-62) 


(9-63) 


3 .   Solution  of  the  Moment  Equations  which  Govern  the 

Growth  of  the  Fluctuations 

The  growth  equations  are  solved  in  Appendix  K. 
The  interested  reader  might  note  the  differences  between 
the  growth  equations  for  the  pressure  free  model  and  the 
radiative  pressure  model.   The  most  significant  difference 
is  that  the  curvature  constant  k  is  in  general  not  negligi- 
ble in  the  region  where  the  matter  density  exceeds  the 
radiation  density.   Therefore,  in  order  to  avoid  difficult 
integration  problems,  the  growth  equations  are  solved  only 
for  the  k  equal  zero  universe.   The  advantage  for  doing 
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this  is  that  the  results  can  be  expressed  in  terms  of  simple 
functions.   (It  seems  that  such  problems  presented  by  a 
nonvanishing  k  are  inherent  to  nonuniform  models.   Sachs 
and  Wolfe  [24],  for  example,  are  also  led  to  consider  only 
Euclidean  models  for  the  same  reason.)   Nevertheless,  the 
results  for  k  equal  to  zero  should  have  some  relevance  to 
k  equal  to  plus  and  minus  one  universes.   This  is  because 
the  dynamics  of  the  universe  are  now  only  slightly  modified 
by  curvature  effects  [25]  and  the  solutions  to  the  growth 
equations  ought  to  therefore  be  approximately  correct  for 
such  models  also. 

With  k  equal  to  zero,  the  fluctuation  free  solutions 

2 

for  Kp  and  R'   given  in  Appendix  E  are 

3a, 


Kp  =  - 


R 


and 


R 


2  _   d 
R 


(9-64) 


where  ad  is  an  integration  constant  different  from  ar  which 
appears  in  the  analogous  equations  for  the  radiation  filled 
universe . 

Then  the  growth  equations  become 


R28rr<t2>  +  |  R8R<x2>  -  3<x2> 


and 


f2 


,% 


R3<O0t)  > 


8r<(30t)2>  +  ^  <(30t)  >  = 


3ad 


—  8r<t2> 
R3   R 


(9-65) 


(9-66) 
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The  solutions  to  these  equations  are  given  in  Appen- 
dix K. 

<(V)2>      -1      -7/2      -6 

ra =  Z0R   +  F0R     +  D0R  (9-67) 

<T2>  =  IqR2  _  2  PqR-1/2  +  |  DqR-3  t  (9_68) 

4 .   Special  Submodels 

Again,  as  in  the  radiation  model,  only  the  following 
three  models  are  considered. 

(1)  <t2>  =  IQR^  ; =  IqR"1  (9-69) 

2 

(2)  <T*>  .  -  I  F0R"1/2  ;  *    "'    *  -  F0R"7/2       (9-70) 

2 

j,  ,         <(80t)    >  c 

(3)  <T2>    =    I  D0IT^     ; =    D0R"b    .  (9-7D 

The  same  difficulty  occurs  with  the  F0-model  as 
with  the  F^^-model  of  the  previous  section.   That  is,  it 
is  impossible  to  find  a  nonzero  F~  with  a  sign  that  is 
compatible  with  both  of  the  Eqs.  (9-70)  (a,  is  always  greater 
than  zero).   Hence,  only  the  In  and  Dn  models  are  considered. 
It  is  to  be  noted  that  the  In-model  is  a  model  for  which 
the  nonunif ormity  increases  with  increasing  R  and  the  DQ- 
model  is  a  model  for  which  it  decreases  with  increasing  R. 

If  Eqs.  (9-69)  and  (9-7D  are  substituted  into 
Eqs.  (9-60)  and  (9-6l)  the  following  equations  are  obtained. 
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(In  -  fluctuations) 


kp  +  57  B(IQ)R 


-1 


3R 


R 


3  B(I  )R 


-1      2R"   R 


R     R' 


(9-72) 


(Dn  -  fluctuations) 


Kp    +    22    B(D    )R 


-6  3R' 


ir 


-    2    B(DQ)R 


-6  2R"        R 


R  R' 


(9-73) 


where 


B(I0)    = 


I0ad 


and 


B(DQ)    = 


Vd 


99-7*0 


be 


The  solutions  of  Eqs .  (9-72)  and  (9-73)  turn  out  to 


(I.  -  fluctuations) 

3a 


Kp  =  - 


d   105  B(V 


(9-75) 


R' 


R 


R'2  =  /"  \  B(Io)R 


(9-76) 


(Dn   -   fluctuations) 

3ad 

Kp     =     -     — 


20    B(DQ) 


R 


R" 
R 


R 

2    B<V     1 
3  Rl 


(9-77) 


(9-78) 
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If  the  same  classification  scheme,  which  is  depicted 
in  Figs.  1  and  2,  is  used  to  classify  these  solutions  that 
is  used  to  classify  the  solutions  for  the  radiation  model, 
then  it  is  easily  found  that  the  I~  -  fluctuations  are  of 
type  K  and  the  D„  -  fluctuations  are  of  type  M.   The  discus- 
sion of  this  model  and  the  radiation  model  is  given  in  the 
next  chapter. 
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X.   SUMMARY  AMD  CONCLUSIONS 

A.   VALIDATION  OF  PROCEDURE 

Before  discussing  the  results  of  this  study  regarding 
the  development  of  a  cosmological  model  to  demonstrate  the 
applicability  of  a  moment  hierarchy  scheme  to  perturbation 
problems  of  general  relativity,  it  would  be  well  to  show 
that  the  solutions  obtained  are  compatible  with  ones  obtained 
by  a  completely  different  procedure. 

Well  known  in  the  literature  is  the  work  by  Sachs  and 
Wolfe  [24]  concerning  the  anisotropy  of  the  microwave  3°K 
background  radiation  due  to  perturbations  of  geometry.   In 
this  study  a  scheme  is  presented  for  analyzing  small,  linear 
perturbations  away  from   an  unperturbed,  k  equal  zero 
Robertson-Walker  metric.   The  method  of  solving  for  the 
perturbations  is  to  first  take  the  spatial  Fourier  transform 
of  the  perturbations,  solve  them  in  k  space,  and  then  trans- 
form them  back  into  position  space.   They  do  this  for  both 
the  dust  and  radiation  models. 

For  the  dust  model  (pressure  free)  they  find  two  kinds 
of  terms.   In  both  types,  6p  decreases  although  6p/p  increases 
for  one  and  decreases  for  the  other  (in  their  notation,  p 
is  the  mass  density  and  6p  is  its  small  perturbation) . 
Furthermore,  the  characteristic  time  for  their  relative 
increasing  solution  is  of  the  same  order  as  the  character- 
istic time  for  the  background  metric.   They  report  their 
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solutions  for  6p  to  be  of  the  form 

_q/2     -2 
6p  =  DR  y/d  +  IR  (10-1) 

where  D  and  I  are  functions  only  of  the  spatial  coordinates. 
Because  p  goes  as  I/R  , 

&■   =  DR"3/2  +  IR  (10-2) 

where  the  constant  of  proportionality  is  absorbed  into  D 
and  I. 

As  shown  by  Eq .  (9-57),  the  ratio  of  6p/p  must  be  propor^ 
tional  to  t.   Hence,  according  to  Eq .  (10-2),  x  must  also 
have  the  form 

x  =  DR~3/2  +  IR.  (10-3) 

Thus,  t2  must  be  given  by 

x2  =  D2R~3  +  (DI)R~1/2  +  I2R2  .  (10-4) 


But  Eq.  (9-68)  says  that 

<t2>  =   I  DQ  R-3  +  —^-   R-1/2  +  IQR2  .        (10-5) 

It  is  seen  that  the  exponents  of  each  term  agree  with 
those  of  Sachs  and  Wolfe.   This  agreement  is  noteworthy 
because  Eqs .  (10-4)  and  (10-5)  are  produced  by  two  very 
different  procedures. 

The  determination  of  the  characteristic  time  TT  for  the 
In  (increasing)  perturbation  term  x  can  be  obtained  from 

2 

the  ratio  of  <x2>  to  <Oqx)  >.   That  is, 
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2      !      <  T2>  ^°"6) 

T    =  —        2 
1     c2  <(30t  )  > 

wherec  is  the  speed  of  light.   Using  Eq.  (9-69) 

(10-7) 


2      R 
T 


1    ***& 


.  «..   ^mp  T  for  the  background  geometry  can 
The  characteristic  time  1R  tor  w 

*.4   ~f  R2  to  (R1  )2  •   That  iSj 
be  found  from  the  ratio  of  R   to  ^ 

.   R2  (10-8) 

2     IK 
m     =  

B    c2  R'2 

.  ^   tr        fa   61H      T   2   becomes 
With  Eq.     (9-^; >    XB 


T 


(10-9) 
R3 


B  c*ad 


■  t    /  nhiPh  is  proportional  to  t  is  pre- 
indeed,  the  ratio  6p/p  which  is  pr  P 

a    i  to  take  place  on  the  same  time  scale 
dieted  by  the  model  to  take  px* 

*r  *v     tn  T  is  one . 
as  the  background  because  the  ratio  of  T,    TB 

•  4  ^-p  soreement  with  Sachs 
Thls  constitutes  a  second  poxnt  of  agreem 

and  W°lfe'  ,    „i„  for  the  dust  models. 

The  above  comparisons  are  made  only  for 

•,on  between  radiation  models  can  be  made  an  a 
Rut  a  comparison  oeu/ecu 

Sachs  and  Wolfe  indicate  that  they  neglect 
similar  manner.   Sachs  an  corres- 

M  in  their  derivation.   Hence,  the  corr 
transport  processes  in  their 

,.-,  of  this  work  is  the  "zero  leakage  model 
ponding  model  of  tnis 

..  h  y    and  6  are  both  equal  to  one. 
WhlC  rtpl   Sachs  and  Wolfe  report  that 

For  their  radiation  model,  Sachs 

(10-10) 

6p  =  DR" 
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so  that 


P 


2  -4 
=  D  R 


(10-11) 


where  D  is  a  function  only  of  the  spatial  coordinates. 

The  general  solution  for  <T2>  given  by  Eq .  (9-36)  with 
X  and  3  equal  to  one  is 


<T  >  = 


¥2r*  + 


D 


1/3  R-H 
Tj-  R 


(10-12) 


The  Pj/o  solution  is  omitted  from  the  general  solution  of 
Eq.  (10-12)  because,  as  stated  in  the  text,  certain  sign 
problems  are  associated  with  it. 
But  Eq.  (9-24)  indicates  that 


<id> 


=  -  2<t2> 


(10-14) 


This  indicates  that 

.  2 


That  is, 


d 
IPJ 


d 
IPJ 


a  t 


'I 

1/3 

R^    + 

^l/sl 

It 

I                       J 

R 


-4 


(10-15) 


(10-16) 


The  comparison  between  Eqs .  (10-11)  and  (10-16)  indicates 
that  although  there  is  disagreement  concerning  the  existence 
of  an   R   term,  there  is  agreement  concerning  the  R   term. 
It  is  not  possible  to  assess  the  reason  for  this  discrepency 
owing  to  the  omission  of  the  solution  procedure  from  Sach's 
and  Wolfe's  paper.   However,  it  is  significant  to  this  work 
that  with  all  the  possibilities  for  the  values  of  the 
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exponents  due  to  the  many  possibilities  for  \   and  3,  there 
is  an  agreement  concerning  the  existence  of  an  R   term. 

Thus,  there  is  a  major  agreement  between  the  functional 
forms  of  both  the  k  equal  to  zero,  dust  and  radiation  filled 
models.   Because  this  agreement  is  accomplished  by  two  very 
different  perturbation  techniques,  support  is  given  to  the 
method  of  this  work  for  treating  perturbations  by  a  hierar- 
chy of  moment  equations. 

It  should  again  be  mentioned  that  the  agreement  of  the 
radiation  model  is  contingent  upon  the  selection  of  the 
"zero  leakage"  values  for  X   and  3  because  Sachs  and  Wolfe 
assume  that  the  velocities  are  zero.   It  is  reasonable  to 
ask  how  valid  is  this  assumption.   In  the  derivation  of  the 
model,  the  attempt  was  made  to  include  the  effects  of  the 
velocities  by  the  introduction  of  the  parameters  X    and  3 
whose  magnitudes  are  estimated  in  Appendix  H  to  be  of 
order  one. 

Figures  3  and  4  give  the  types  of  solutions  that  can  be 
expected  for  all  possible  combinations  of  X    and  3 .   In  the 
case  of  the  In  /r.   -   type  perturbation,  the  model  classifica- 
tion for  X    and  3  values  of  order  one  is  clearly  type  K.   But 
for  the  D,/:.  -  type,  both  the  K  and  M  -  type  solutions  are 
possible.   In  fact,  the  zero  leakage  values  of  X   and  3  are 
observed  to  lie  on  the  boundary  between  the  K  and  M  -  type 
solutions.   Thus,  if  a  D^-  -  type  model  is  being  considered, 
the  neglect  of  the  velocities  is  not  trivial  for  very  differ- 
ent dynamical  situations  could  occur  even  if  the  velocities 
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are  extremely  small.   But  for  an  I  .   -  type,  the  neglect  of 
the  velocities  does  not  change  the  solution  type.   In  this 
case  the  approximation  of  zero  velocities  seems  to  be  justi- 
fied.  Indeed,  Sachs  and  Wolfe  consider  only  this  type. 

However,  X  and  3  not  only  affect  the  classification  but 
also  determine  the  exponent  (v/hich  is  j  +  4  for  D,  ,~  and  i  +  4 
for  L  ,J  of  R  in  the  expression  for  <t2>.   If  Figure  5 
the  lines  of  constant  i  and  j  are  given  in  X-3  space.   Prom 
these  curves,  it  appears  that  for  X  and  3  of  order  one, 

-3  5  j  +  4  <  -  8 

3  S  i  +  4  <  8  .  (10-17) 

In  the  case  of  X  and  3  equal  to  exactly  one,  j+4  is  equal  to 
-4  and  i+4  equal  to  4.   These  estimates  are  conservative. 
Perhaps  a  more  reasonable  estimate  would  be  to  assume  that 
the  likelihood  of  X  and  3  both  being  on  the  large  side  of 
their  order  one  magnitude  is  small.   Then  the  upper  limit 
might  be  pushed  down  to  perhaps  -6  and  6  respectively. 

Thus  it  seems  that  the  zero  velocity  assumption  can 
reasonably  be  questioned  on  the  grounds  that  the  exponent  of 
<t2>  could  be  different  enough  from  the  zero  velocity  expo- 
nent that  a  many  order  of  magnitude  change  of  R  could  produce 
a  value  of  <t2>  that  is  much  larger  (or  smaller)  than  the 
corresponding  zero  velocity  value.   This  can  be  seen  more 
quantitatively  as  follows.   Consider,  for  example,  a  radia- 
tion model  with  an  I1/^  -  type  perturbation.   At  some  value 
of  R  equal  to  R, 
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1± 

<iZ>l  =  ~W~  Rl+i4  (10-18) 

where  the  sub  and  superscript  1  refers  to  the  i-value  of  the 
exponent.   If  it  is  assumed  that  the  ratio  <t2>./<t2>  ±s 
one  at  R,  then  this  ratio  is  in  general  different  at  other 
values  of  R.   That  is,  this  ratio  is  a  measure  of  the  effect 
of  a  given  combination  of  A  and  3  upon  the  nonunif ormity 
in  terms  of  the  "zero  velocity"  nonunif ormity .   This  ratio 
is  easily  found  to  be,  for  an  arbitrary  R 


'R1 


-2>o 


R 


(10-19) 


where  i  =  -  4  +  2/2(A+3)  . 

Thus,  unless  A  and  3  are  such  that  i  is  near  enough  to 
zero  so  that  <t2>  is  of  order  <t2>  ,  the  effects  of  the 
expansion  through  many  orders  of  magnitude  upon  the  nonuni- 
formity  compared  to  the  zero  velocity  nonunif ormity  is  signif- 
icantly increased  as  seen  from  Eq.  (10-19).   If,  however, 
R  changes  at  most  by  only  a  few  orders  of  magnitude  in  a 
given  problem,  the  discrepency  is  not  as  serious  and  may  be 
tolerable  if  i  is  close  enough  to  0.   Hence,  although  the 
zero  velocity  approximation  is  extremely  convenient  mathe- 
matically, its  validity  does  not  appear  to  be  on  solid  footing 

B.   DESCRIPTION  AND  COMPARISON  OF  MODELS 

As  stated  in  Section  A  of  this  chapter,  there  are  two 
intrinsically  different  types  of  nonunif ormities  associated 
with  both  the  dust  and  radiation  filled  models.   One  is  an 
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increasing  (I-type)  and  the  other  is  a  decreasing  (D-type) 

2 

nonuniformity.   If  the  ratio  of  <t2>  to  <(8qt)  >  is  computed 
for  all  four  models  (two  for  each  a  value)  in  the  same  way 
as  done  previously  so  that  a  characteristic  time  might  be 
estimated  for  the  relative  perturbations,  then  it  is  found 
that,  for  all  models,  this  characteristic  time  is  of  the 
same  order  as  the  characteristic  time  for  the  average 
background . 

Because  of  the  existence  of  the  two  D  and  I-type  pertur- 
bations, it  is  natural  to  ask  which  kind  characterizes  our 
universe.   This  question  can  only  be  answered  through  obser- 
vation.  Although  such  observations  are  difficult  to  perform, 
there  seems  to  be  some  suggestion  that  the  I-type  of  pertur- 
bation is  the  one  particular  to  our  own  universe.   One 
reason  why  this  might  be  true  is  that  the  average  background 
metric  might  be  expanding  more  rapidly  than  the  superc luster 
to  which  our  galaxy  belongs.   A  second  reason  is  that,  in  the 
case  of  D-type  perturbations,  galaxies  might  be  discarded 
from  their  otherwise  clustered  configurations  \_2h~\  . 

Granting  the  I-type  dust  model,  it  is  possible  to  estimate 
the  parameter  In  which  governs  the  amplitude  of  the  metric 
perturbation  <t2>.   This  can  be  accomplished  by  using  the 
data  reported  by  Partridge  and  Wilkinson  [19]  on  the  aniso- 
tropy  of  the  so-called  cosmic  blackbody  radiation.   They 
report  that  the  ratio  6T/T  is  of  the  order  10   where  T  is 
the  temperature  of  the  radiation  (3°K)  and  6T  is  the  fluctua- 
tion of  T  with  angle  of  observation.   Sachs  and  V/olfe  [24] 
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give  a  very  simple  expression  that  relates  6T/T  to  <SP  /p 
where  p0  is  the  current  mass  density. 


T    2   p0  CHRL) 


(10-20) 


where  H^  = 


V 


R   R, 


('  denotes  derivative  wrt .  x°),  R   is  the 


current  value  of  R,  and  L  is  the  characteristic  length  of  the 
perturbations.   By  Eq.  (9-57),  Eq.  (10-20)  can  be  written  as 


_2   .  16 
9 


(HRL)* 


6T 


(10-21) 


The  constant  I   can  be  estimated  by  means  of  Eqs.  (9-69) 
and  (10-21) .   That  is, 


R, 


ill  f  ii 

T  I  I  9 


^HRR0)4 


(10-22) 


where 


R 


(10-23) 


If  HRR   is  taken  to  be  one  and  if  the  factor  16/9  is  neglect- 
n  0 

ed  since  it  is  of  order  one,  Eq.  (10-22)  becomes  simply, 


f6T 


0 


R 


^2 


rr 


(10-24) 


If  the  fluctuations  are  assumed  to  be  the  size  of  about 
1000  MPC  [24,  31]  (one  MPC  is  equal  to  3.26  light  years) 
so  that  e  is  of  order  ten  (R   ~  1029  cm),  then  for  Partridge 


and  Wilkinson's  estimate  of  6T/T,  Eq.  (10-24)  gives  for  I 


0 
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I0  -  10~58cm  2    .  (10-25) 

Using  this  value  for  IQ,  it  is  interesting  to  determine 

the  effect  a  term  of  this  strength  would  have  upon  the  average 

dynamics  if  it  is  assumed  that  the  long  wavelengths  are 

important  to  the  expansion  characteristics  [35].   It  is  well 

known  that  the  k  =  0,  unperturbed  uniform  universe  model 

expands  forever  [2].   (This  can  be  seen  easily  by  integrating 

Eq.  (9-64)  to  find  R  a  t2/3.)   But  owing  to  the  additional 

source  term,  the  dynamics  of  the  model  is  perturbed.   As 

pointed  out  earlier,  the  solution  for  the  k  =  0,  I-.  model 

is  type  K.   Because  this  means  that  R  is  perturbed  below 

the  uniform  solution,  (see  Fig.  1),  the  question  naturally 

arises  whether  or  not  there  exists  a  maximum  value  for  R. 

To  answer  this,  R'   must  be  set  to  zero  in  Eq .  (9-76)  and  a 

positive  value  of  R  greater  than  R0  must  be  determined. 

this  is  done,  such  a  value  R    is  found  such  that 
'  max 


R 

max 


'16  U/2 


3J0 


(10-26) 


Using  Eq.  (10-25) , 

R     »  1029  cm  .  (10-27) 

max 

That  is,  if  the  observed  anisotropy  of  the  blackbody  radia- 
tion is  due  to  perturbations  of  space  as  large  as  super- 
cluster  distances  and  are  of  the  IQ  -  type,  then  the 
"effective  energy"  of  these  longwave  fluctuations  is 
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sufficient  to  prevent  an  unlimited  expansion  of  a  k  =  0, 
dust  filled  universe  which  would  otherwise  expand  forever. 
It  should  be  pointed  out  that  this  model  is  not  valid 
for  all  values  of  R.   For  small  R,  the  radiation  density, 

which  goes  as  R   ,  overwhelms  the  mass  density,  which  goes 

_3 

as  R   .   At  this  point,  it  is  no  longer  valid  to  speak  of  a 

dust  model  and  it  is  proper  to  use  the  radiation  model  for 
smaller  R-values  [15].   On  the  other  hand,  for  large  R, 
<t2>  approaches  one.   As  it  does  so,  the  assumption  of  small 
amplitude  fluctuations  which  was  used  to  truncate  the  moment 
hierarchy,  becomes  invalid.   The  point  at  which  this  happens 
can  be  easily  determined.   With  <x2>  equal  to  one,  Eq.  (9-69) 
gives 

R   =  ^-7^  (10-28) 

c    (   }l/2 

where  R   is  that  value  of  R  for  which  <t2>  is  equal  to  one. 
c 

If  Eq.  (10-28)  is  compared  to  Eq.  (10-26),  it  is  found 

that  R  <R    which  indicates  that  at  the  "turning  point," 
c   max 

the  fluctuations  are  as  large  or  larger  than  the  average 
metric.   However,  it  should  be  noted  that  because  the  model 
is  not  applicable  at  the  turning  point  itself,  it  may  not 
therefore  be  justifiably  used  to  determine  RmQ.   This  diffi- 
culty  is  not  serious  though  since  the  ratio  of  the  "extrapo- 
lated value"  of  R  ov  given  by  Eq.  (10-26)  to  R   is 

max  i* 


rcax  =  H  (10-29) 

Rc    /3 


90 


which  is  of  order  one.   Thus  the  real  value  of  R    should 

max 

still  be  of  the  order  of  1029  cm  as  given  by  Eq .  (10-27) 

since  the  theory  breaks  down  close  to  the  estimated  value 

of  R    .It  should  finally  be  noted  that  if  other  models 
max  J 

are  considered  which  also  have  "turning  points,"  (either  as 
a  minimum  or  maximum)  the  weak  fluctuation  assumption  breaks 
down  just  before  the  turning  point  occurs.   Thus,  if  the 
physics  of  an  assumed  turning  is  to  be  accurately  discussed, 
a  strong  fluctuation  model  should  be  used. 
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APPENDIX  A 
The  Reynold  Stresses 

The  Navier-Stokes  equation  for  an  incompressible  fluid 
is  [12] 

±I+V.VV=--Vp+  vV2V  (A-l) 

a  t  p 

where  V  and  P  are  the  velocity  and  pressure  fields  relative 
to  a  fixed  coordinate  system  and  v  is  the  ratio  n/p  where 
n  is  the  viscosity  coefficient  and  p  is  the  mass  density. 
For  incompressible  flow,  p  is  a  constant  throughout  all 
space . 

In  order  to  gain  some  insight  into  the  perturbation 
technique  that  is  presented  in  this  paper,  the  reader  should 
look  upon  the  Navier-Stokes  equation  as  if  It  were  the 
Einstein  Field  equations.   The  left  hand  side  of  Eq .  (A-l) 
is  analogous  to  the  geometrical  side  of  the  field  equations 
and  the  right  hand  side  is  analogous  to  the  mass  side 
(that  is,  the  side  which  contains  the  energy  momentum  ten- 
sor).   The  velocity  field  V  plays  the  same  role  in  the 
Navier-Stokes  equation  as  does  the  metric  tensor  g    in  the 
Einstein  field  equations. 

In  fluid  turbulence  theory,  the  following  decompositions 
are  assumed  [4].   (In  this  appendix-,  the  equation  after  the 
semicolan  refers  to  the  analogous  equation  in  the  text.) 

V  =  <V>  +  v  (A-2);(4-9) 
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P  =  <P>  +  p  (A-3);(i]_9) 

where  <  >  denotes  "average." 

If  Eqs.  (A-2)  and  (A-3)  are  substituted  into  Eq .  (A-l), 
the  following  equation  results 

2<|>  +  ||  +  <^>  .   v<v>  +  *  .   v<V> 

+  <V>  •  Vv  +  v  •  Vv 

=  ,  ±   y<P>  -  -  Vp  +  vV2<V>  +  vV2v.         (A-4) 
P         9 

In  fluid  turbulence  theory  it  is  assumed  [4]  that  if 
A  is  a  random  field, 

<|f>=|T<A>  (A-5);CB-2) 

and 

<9A_o    =    L_  <A>     .  (A-6);(B-2) 

8x  dx 

Because    of   Eqs.     (A-2)    and    (A-3), 

<v>   =   0  (A-7);(4-l4) 

and 

<p>    =    0    .  (A-8) 

Then,   Eq.  (A-4)  becomes 

-> 

?£XL  +  <V>  •  V<V>  =  -  -  V<P> 
9 1  p 

+  vV2<V>  -  <v  •  Vv>  .         (A-9);(6-8),(6-9) 

It  should  be  noted  that  Eq.  (A-9)  is  symbolically  iden- 
tical  with  Eq.  (A-l)  with  V  and  P  replaced  with  <V>  and  < P> . 
The  only  difference,  however,  is  the  appearance  of  the  extra 
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source  term  from  the  "velocity  side"  which  results  from  the 
nonlinear  V  •  VV  term.   It  is,  in  general,  nonvanishing. 

This  extra  term  is  called  the  Reynold  Stress.   It  can 
be  interpreted  as  an  additional  source  term  (in  addition 
to  the  pressure  and  viscous  sources)  which  acts  on  the 
average  flow  field  (<V>).   Such  an  interpretation  can  be 
made  because  the  left  hand  side  of  Eq .  (A-9)  is  symbolically 
identical  with  the  left  hand  side  of  Eq.  (A-l). 

In  order  to  obtain  the  functional  form  of  the  Reynold 
Stresses,  a  hierarchy  of  equations  can  be  built  up  from  the 
Navier-Stokes  equation  of  Eq .  (A-l)  by  multiplying  it  with 
various  orders  of  V  and  averaging.   But  the  nonlinear  term, 
as  already  seen,  always  produces  a  moment  of  one  order 
greater  than  the  order  of  the  moments  which  comes  from  the 
other  linear  terms.   Thus, there  is  always  one   more  unknown 
than  there  are  necessary  equations  available  to  solve  the 
system  of  equations.   This  problem  is  the  well  known  closure 
problem  of  turbulence  theory  which  is  discussed  in  the 
text  [3J. 
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APPENDIX  B 
Definition  of  the  Averaging  Operation 

There  are  various  definitions  of  average  which  are  pos- 
sible.  Many  theorists  who  work  in  the  field  of  fluid  turbu- 
lence use  an  integral  expression  for  their  definition  of 
average.   For  example,  Curie  and  Davies  [4]  use  a  time 
integral.   This  is  an  acceptable  definition  for  experimental 
work  but  as  pointed  out  by  Lumley  and  Panofsky  [13]  }  such 
an  integral  expression  is  not  always  easy  to  define  in  a 
meaningful  way.   They  suggest  that  an  ensemble  average  is 
probably  the  most  basic.   Thus,  the  notion  of  an  ensemble 
average  is  taken  as  the  definition  for  the  averaging  process 
<  >  for  this  model. 

Suppose  D(x^,£)  is  a  stochastic  process  where  x^  are  the 
four  space-time  coordinates  and  £  is  the  ensemble  parameter 
which  designates  D  as  a  possible  outcome  for  an  experiment 
(see  also  Ch.  IV).   Then  <D>  is  defined  as 

D(xy;  5.)  +  ...  +  D(xy;  £  ) 
<D>  =  lim —  (B-l) 

where  <D>  is  only  a  function  of  xy  and  is  not  dependent  on 
£.   This  is  the  same  definition  as  that  used  by  Lumley  and 

Panofsky. 

Because  this  definition  involves  a  summation,  the  follow- 
ing equation  is  assumed  valid. 
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d    <D>  =  O  D>   .  (B-2) 

Suppose  it  is  desired  to  know  the  probability  of  finding 
D(x  )  falling  below  a  certain  value  6.   This  can  be  deter- 
mined by  counting  the  proportion  of  members  of  the  ensemble 
for  which  D(x  )  1  6.   From  this  a  function  F  can  be  defined 
as  follows 

F(6;xy)  =  P{D(xy)  1    6}  (B-3) 

where  P  stands  for  "Probability  that."   The  probability 


density  function  is  defined  by 

8F(6;xP) 
f(D;xy)  =  r6 .  (B-4) 

Then  the  average  of  any  function  h[D(x  )]  can  be  written  as 

oo 

<h[D(xy)]>  =  /  h(6)f(6,xy)d6  .  (B-5) 
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APPENDIX  C 
An  Implication  of  Spatial  Stationarity  of  the  Random  Field 

An  implication  of  statistical  spatial  stationarity  is 
rather  easy  to  state  because  of  the  assumed  form  of  the  line 
element.   On  a  hypersurface  of  constant  x° ,  dx°  is  zero  so 
that  by  Eqs.  (4-5)  and  (4-6) 


ds   o 


V(x\xJ';0 


dc 


(C-l) 


1  + 


k  r' 
IT 


v, 


It  is  pointed  out  in  Chapter  IV  that  RM  has  the  meaning  of 
"local  radius  of  curvature."   It  is  also  a  stochastic  process 
which  means  that  at  each  point  x   of  space  time  it  has  a 
probability  density  function  associated  with  it  (see  Appen- 
dix B)  . 

In  order  to  bring  about  a  simplification  of  the  proposed 
model,  consider  the  following  assumption.   Although  the 
hyperspaces  of  constant  x°  are  clearly  nonuniform,  and  under 
go  various  expansions  at  each  point,  it  is  assumed  that  the 
probability  of  finding  a  given  expansion  at  any  point  is 
independent  of  position  on  that  hyperspace.   This  require- 
ment is  then  one  of  statistical  spatial  stationarity,  a 
concept  which  is  familiar  in  fluid  turbulence  theory.   This 
implies  that  the  probsbility  of  measuring  a  given  value 

\   2 


R  2  between  say  R  2  and  R  2  +  dR  *  ought  to  be  independent 


M 


!•: 
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of  the  coordinates  on  this  spatial  hypersurface  of  constant 
time  x°  because  this  random  function  contains  all  the  expan- 
sion characteristics.   Equivalently ,  the  probability  of 
finding  t  between  T  and  T  +  dT  ought  to  also  be  independent 
of  the  spatial  coordinates.   Hence,  the  probability  density 
function  f(T;xy)  defined  by  Eq.  (B-4)  must  be  such  that, 

f  (T;x°,xJ  )  =  f  (T;x°,xJ*  ') 

i     i  ' 
for  every  xJandxJ   where  j  =  1,  2,  3.   This  can  only  mean 

that  [18] 

f(T;xy)  =  f(T;x°)  (C-2) 

That  is,  f  can  only  be  a  function  of  the  coordinate  x°  and 
t.   But  Eq.  (B-5)  says  that  for  any  function  h  of  the  process 
x(xy;£)> 

00 

<h(T)>  =  /  h(T)  f(T;x°)dT  .  (C-3) 

—00 

Hence,  <h(x)>  can  only  depend  upon  x°  and  not  upon  r,  6,  <j> . 
It  follows  that 

8.<h(T)>  =  0  (C-4) 

J 

where  j  =  r,  6,  cj) . 
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APPENDIX  D 
The  Form  of  the  Energy-Momentum  Tensor 

The  general  form  of  the  energy-momentum  tensor  for  a  per- 
fect fluid  is  given  by  [1,  29] 

T^v  =  DU^UV  +  £5-  (UUUV  -  gyV)  (D-l) 

c 

where  U^  =  dxy/ds  and  D  and  P  are  the  proper  hydrostatic 
density  and  pressure.   That  is,  these  are  the  density  and 
pressure  which  would  be  measured  by  an  observer  moving 
instantaneously  with  the  fluid.   Tolman  [29]  defines  a  per- 
fect fluid  to  be  a  mechanical  medium  incapable  of  exerting 
transverse  stresses  in  the  proper  frame.   Prom  this  defini- 
tion follows  Eq.  (D-l). 

Tolman  [28]  shows  that  the  energy-momentum  tensor  for  a 
disordered  electromagnetic  field,  which  is  expressed  in  terms 
of  the  electric  and  magnetic  field  strengths,  can  be  viewed 
as  a  perfect  fluid  for  which 

F  -  j  c2D.  (D-2) 

However,  if  the  pressure  is  due  to  the  random  motion  of 
particles  of  cold  dust,  then  the  equation  of  state,  which 
is  well-known  to  kinetic  theory,  is  given  by  [1] 


P  =  i^-  p  (D-3) 

3  c2 
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where  V   is  the  root-mean-square  random  velocity  of  the  dust 
particles.    But  this  velocity  is  much  less  than  c  so  the 
pressure,  for  all  practical  purposes,  is  essentially  zero. 
The  two  models,  therefore,  can  be  written  together  in 
one  equation. 


P  =  ac2D 


CD-4) 


where    a  =  1/3  if  universe  is  radiation  dominated 
a  =  0    if  universe  is  matter  dominated  . 


Now,  the  expression  for  the  mixed  energy  momentum  tensor 
T^   is  needed.   This  is  easily  found  by  means  of  the  metric 


tensor  g 


yv 


Ty   =  DUyU 
v 


+  —  (UyU   - 


6M  ) 
v 


(D-5) 


where 


U   =  g   U 


a 


(D-6) 


In  order  to  evaluate  Eq.  (D-5),  U   and  U   need  to  be 
determined.   This  can  be  done  as  follows.   Consider  a  fluid 
element  which  moves  along  the  trajectory  xJ(x  ).   Along 
such  a  trajectory, 


ds2  =  (dx°) 


1  -  e 


M  da 


(dx°) 


(D-7) 


With  the  definition 


u  =  e 


M/2 


dx 


(D-8) 
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equation  (D-7)  becomes 


ds 
dx° 


+  (1  -  u2) 


1/2 


(D-9) 


Notice  that  if  the  motion  of  the  fluid  element  is  along  a 
null  geodesic,  u2  must  be  equal  to  one.   In  order  to  keep 
the  interval  of  such  trajectories  timelike  (interval  real) 
so  that  the  local  velocity  of  the  fluid  element  does  not 
exceed  the  light  velocity,  it  is  obvious  that  from  Eq .  (D-9) 


0  <  u2  <  1  . 


(D-10) 


It  is  natural  to  identify  (uc)  with  the  local  velocity  of 

M/2 
the  fluid  element  at  the  point  x^  because  e   da  is  the 

physical  distance  that  the  fluid  element  moves  in  a  time 

—  dx° .   Now  define 
c 


Y 


(1  -  u2) 


1/2 


(D-ll) 


so  that 


v  o 


u° 

ds 

- 

Y  • 

Therefore , 

uJ 

dxJ" 

_ 

dxJ*  dx° 

dxJ 

ds 

dx°  ds 

dx° 

Y 


(D-12) 


(D-13) 


But  Eq.   (D-8)  can  be  written  as 


u2  =  e 


M 


dr 


dx1 


+  r- 


de 


dx 


+   r2sin2 


r  ^    ^  2 


dx° 


.   (D-14) 


If  the  quantities  ur,  uQ  and  u   are  defined  such  that 
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M/2   dr 
u„   =   e 


dx 


M/2  d9 
u     =  re  — 


e  dx° 

u      =   r    sine    eM/2   ^L_  (D_15) 

*  dx° 

then 

u2    =   u    2    +   u   2    +   u    2  (D-16) 

r  6  <|> 

It  is  natural  to  identify  u.  as  the  j     component  of  the 
local  velocity  at  x  .   Hence, 

U°    =   y 

TIi         v    -M/2 
U1    =    Ye  u 

-M/2 
o               e 
U2    =    Y un 


-M/2 
U3    =   Y  .    nu±    .  (D-17) 


r    sin6     <$> 
By   Eq.     (D-6), 


u 

0 

— 

Y 

U 
1 

= 

- 

v   M/2 
Ye        u 
r 

u 

2 

= 

- 

v      M/2 
Yre        u0 

U      =   -   Yr    sine    eM/2uA    .  (D-18) 

3  <t> 
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Then  Eq.  (D-5)  becomes, 


T°      =   D   +  —  u2 

0  2 

c  z 


Ti      =    -   Du  2-   —   (1    +   u  2-   u2) 


T2      =    -   Du  2-   L.   (1+   u  *'_   u2) 

e      C2  e 


T3      =    -   Du  2-   —   (1    +   u  2-   U2) 

3  4>         2  4> 


*•      -   -   eV„    -  -    (D  +  2-)    eM/2u 

c2  J 


T»      =   -  r2eHT2o    -   -    (D   +  E-)    reM/2u, 

2  0  2  ' 

C 


»  ?  5         M    3  ,  P    v  M/2 

a    =   -   r2sin26    e   T  ft     =   -    (D  +  ~)   r    sine    e   '    u. 

3  0  2  (|> 


Tl,    =   rV,     =   -    (D  +  £j-)    ruruQ 


T  ,    =   r2sin29    t\     =   -    (D   +  — )    r    sine    u   \ix 
3  J  c2  r   <J> 

T23     =    sin26    T32     =    -    (D    +   — )    sine    u   u  (D-19) 

c2  ^ 

where        D    =   D~Y2 

P   =    PY2    . 

To  produce  Eqs.(D-19),  the  following  identities  are  needed. 

Y2  -  1  =  Y2u2 

Y2u?   +  1  =  Y2(l  +  u.2-  u2)   .  (D-20) 
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APPENDIX  E 


m 


Solution  of  Cosmological  Equations  with  a  1/R   Source  Term 


The  two  usual  cosmological  equations  of  general  rela- 
tivity with  zero  cosmological  constant  but  with  a  1/R 
source  term  are  [1] 


Kp  +  Kp 


V 


3ii  _  3RJ. 
R2    R2 


(E-l) 


and 


kp   kPv 


c2    c2 


k_   2RJ_'  _  R^_ 
R2     R    R2 


(E-2) 


where 


V_ 
v  '  Rm 


(E-3) 


P  -  ac2 


au-p 


P   =  a  c^p 
v    v   Kv 


(E-4) 

(E-5) 


K  =  - 


8ttG 


(E-6) 


dx 


The  difference  between  Eqs .  (E-l)  and  (E-2)  is 


(E-7) 


k(p  +   Py)  +^(P+Pv)  +||=  2 

C  A 


'R'' 


R 


(E-8) 


The  derivative  of  Eq.  (E-l)  is 


K(p«  +  p')  =  — —  -  6 

v      R3 


fR'l 

[r  J 

[R'l 

[r  J 

(E-9) 
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The  substitution  of  Eq.  (E-8)  into  Eq .  (E-9)  gives 

(p'+pv'  )  +  ^-   (p  +  pv)  +  3RL  _JLz_  =  o  .      (E-10) 

c2 

If  Eq.  (E-8)  is  multiplied  by  R3,  it  becomes 

t    (P+P  )      , 
((P  +  PV)R3)   +  ~-  (R3)   -  0   .  (E-ll) 


This  is  the  conservation  of  mass  eouation  where  p   and  P 

Hv      v 

are  the  mass  density  and  pressure  of  the  source  respectively 
It  is  valid  for  any  equation  of  state  P(p)  and  P  (p„). 

But  with  Eq.  (E-5),  P   =  a  c2p   and  Eq.  (E-4),  P  =  ac2p 
Eq .  (E-ll)  becomes 

(pR3)'.  +  ap(R3)'  =  -  (PVR3)'  -  avPv(R3)'  .     (E-12) 

The  homogeneous  equation  associated  with  this  equation  is 

(pR3)'  +  ap(R3)'  =  0,  (E-13) 


the  solution  of  which  is 

p  .     A  (E_l4) 

R3(cx  +  1) 

In  order  to  solve  Eq.  (E-12),  Eq.  (E-l4)  may  be  substi- 
tuted into  it  assuming  A  is  now  a  function  of  x° .   Eq.  (E-12) 
becomes 

A'  =  -  Pv'R3a+3  -  3R3a+2R'PvU+txv)   •         CE-15) 

With  p   =  —  ,  Eq.  (E-15)  becomes 
Kv   Tjin  ' 
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But  Eq.  (E-16)  is  easily  integrable.   It  turns  out  that 


A  =  -  CVR3(a+1)-m  +  B 


where 


C  = 


m-3(l+av) 
m-3(l+a) 


(E-17) 


(E-18) 


Finally,  the  substitution  of  Eq.  (E-17)  into  (E-HO  gives 


B      CV 
p   R3(a+1)    Rm 


(E-19) 


The  substitution  of  Eq.  (E-19)  into  Eq .  (E-l)  gives  an 


expression  for  R'  . 


R 


B   +  V(l-C) 


3a+l 


R 


R 


m-2 


-  k 


Thus, 


Kp 


3a 


R3(a+D    Rm 


and 


R 


+   b 


where 


a  =  - 


R3a+1   Rm-2 


KB 
3 


-  k 


f  =  KCV 


b  = 


kV(C-I) 
.  3 


(E-20) 
(E-21) 
(E-22) 

(E-23) 
(E-24) 

(E-25) 


If  V  =  0  (no  source  term),  then  the  solutions  become 
the  usual  cosmological  solutions  for  p  and  R.  For  a  radiation 
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filled  universe  (a  =  1/3)  [29]  , 

3a 

Kp  = (E-26) 

R* 

and 

a 
R'   =  —  -  k  .  (E-27) 

R2 

For  a  pressure  free  universe  (a  =  0),  the  zero  source 
term  solutions  are 

kp  = ±  (E-28) 

R3 


and 


a, 
R'   -  gS.  -  k  .  (E-29) 


It  should  be  noted  that  the  meaning  of  the  constant  a 
is  different  in  the  case  of  the  radiation  model  (denoted 
by  a  )  than  in  the  case  of  the  pressure  free  model  (denoted 
by  aj. 
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APPENDIX  P 

Elimination  of  the  Spatial  Derivatives 

from  the  Dynamic  Equations  for  a  Radiation  Pilled  Universe 

The  terms  involving  the  spatial  partial  derivatives  of 
t  can  be  removed  from  Eqs .  (9-1)  and  (9-2)  by  the  substitu- 
tion of  Eq.  (9-3).   The  result  is 

Kp  '  2  R~  8o<t2>  +  i   <^ot)  >  +  —  <t2> 

R 

=  12k  +  mn  +  i2Rjl  (f_1} 

R2  R  R2 

-    KP    -   |  ~    80<x2>    -    |    32<t2>    +    |   <(80t)2>    +   ~2    <x2> 

2 

4k         3R "        ilR » 
=    —  +   ^-  +   — .  (F-2) 

R2  R  R2 

But  the  right  hand  side  of  Eqs.  (F-l)  and  (F-2)  are  not 
the  expressions  that  would  result  from  substitution  of  the 
average  metric  into  the  Einstein  field  tensor.   In  order  to 
bring  these  equations  into  such  a  form,  Eq .  (F-2)  must  first 
be  multiplied  by  minus  five  and  added  to  Eq.  (F-l).   If  this 
sum  is  multiplied  by  three  eighths,  then  the  desired  form 
results . 

KP  +  f|  8^<T2>  +  I    gi  80<T2>  -  I    <(90T)% 

3  k   _;>   _   3k   3R' 


<T2>  m   _  2£  _  2£1_  .  (F_3) 

R2  R2     R2 
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To  obtain  the  other  equation,  Eq .  (F-l)  must  be  multiplied 
by  minus  5/3  and  Eq.  (F-2)  by  3.   If  these  two  equations 
are  added  and  divided  by  eight,  the  resulting  equation  is 
formed . 

i  £_  <Tt>  .  _  S_  _  Si*  _  £L.     .  cp-io 

R2  R2  R  R2 

2 

The  term  dQ    <t2>  can  be  replaced  by  using  Eq .  (9-*0-   The 
final  pair  of  equations  are 

KP     +     |     K<Td>     +     |    <O0T)2>     +     f    jp     30<T2> 

_   3L<l2>   =  .   3k  _   3R^  (F.5) 

*   R2  R2  R2 

_  <£.  J  K<Td>  _  1  <O0t)2>  +  $£-  s0<x2> 


1    k_  <t2>         _   L_  _    2R^'  _   RJ_ 
^  R2  R2  R  R2 


2 


(F-6) 
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APPENDIX  G 
Elimination  of  the  Spatial  Derivatives 
from  the  Dynamic  Equations  for  a  Dust  Pilled  Universe 

The  elimination  of  the  terms  which  contain  the  spatial 
partial  derivatives  in  Eqs.  (9-48),  (9-^9)  and  (9-50) 
gives 

-  |  kP  -  |  |^-  80<t2>  +  |-  <O0t)  >  -  |  K<Td> 

+  3k<l2>  .  Igk  .  12R^+  15RL'  (G_X) 

R2  R2      R2       R 

-  f  *p  -  IF  vt2>  -  iK<^> +  Hv>2> 

.  5  K<Td>  +  k_  <t2>  =  4k  +  ig/  +  3RJ_'  m       (G_2) 

5  R2         R2     R2      R 

Now,  if  Eq .  (G-2)  is  multiplied  by  minus  five  and  added  to 
Eq .  (G-l)  and  the  resulting  equation  multiplied  by  three 
eighths,  the  following  equation  is  obtained. 

+|K<xd>  -  ^<x2>  =  _^_  3Rll  .     (G.3) 

0  M  r2  R2       R2 

If  Eq.  (G-l)  is  multiplied  by  minus  five  thirds  and 
Eq.  (G-2)  by  three,  and  if  these  two  resulting  equations 
are  added  and  divided  by  eight,  the  following  equation  is 
obtained. 
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1       R'  0  ^00  l,v2  <5 

?  r-  VT  >  ~  TE   so<T  >  +  F  <(  V}  >  +  2T  K<Td> 


1  k 
R 


T, 


<T">  =  - 


k_ 
R2 


2R"   R 


(0-4) 


R    R; 


Now,  the  term  82<i2>that  appears  in  Eqs .  (G-3)  and  (G-4) 
can  be  eliminated  by  using  Eq.  (9-51).   The  result  is 

kP  +  \   <O0t)  >  +  |  |^-  d0<T2>   +   |  «xd> 


R2 


3k  _  3RJ_ 
R2    R2 


(G-5) 


k   <(V)2>  +  ff  ¥-   8o<t2>  +  IT  K<Td> 


1  k    ,      2X   _ 

R2 


k_ 
R2 


R     R2 


(G-6) 
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APPENDIX  H 

The  Dynamics  of  the  Energy  Momentum  Tensor 

for  Disordered  Radiation 

The  dynamics  of  the  disordered  radiation  fluid  is 
described  by  Eq .  (9-10) 

T^v;  v  =  0 

where 

Tyv        Tyv  +   v  Tyo  +  rW  Tav  (H_1} 

v       va       av 

The  components  of  T  *\  are  given  in  Appendix  D,  and  the  com- 


ponents of  r     are  given  by  Eq .  (5-10).   Because  of  the 
a 


vy 

summation  which  appears  in  Eq.  (H-l)  with  respect  to  the 
indices  v,  Eq.  (H-l)  represents  a  set  of  four  differential 
equations.   Adler  [1]  points  out  that  the  equation  with  y 
equal  to  zero  is  the  energy-conservation  equation  while 
the  remaining  three  are  the  relativistic  equivalent  of  New- 
ton's law  of  dynamics  applied  to  an  ideal  fluid.   The  con- 
servation of  energy  is  therefore  expressed  by 

3  Tov  +  pv  Toa  +  ro  Tav  _  Q   <  (H-2) 

v       va       av 

With  the  appropriate  substitutions  from  Appendix  D  and 
Eq.  (5-10),  Eq.  (H-2)  becomes 

3  (D  +  L_  U2)  +  1  (3  M)[(D  +  —  )(3  +  u2)] 
0     c2       2    °        c2 


+  e"M/28„[(D  +  L-  )Uyi  ]  (H-3) 
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r      c2    r  cont. 


+  e  M/2(3  M  +  |)[(D  +  ^)u  ]  cont. 

P     r        c2   r  (H-3) 

+  e~M/2  I  3  [(D  ♦  ?-)uel 

c 

t  e-M/2(  i  3eM  +  £2i_e)[(D  +  Ej)  u9] 

+  e"M/2  r^IFT  V(D  +  h   V 

c 

+  e~M/2  ( i 3-  3.M)[(D  +  — )  uj  =  0  . 

vr  sin  0  <b  2        h> 

c 

If  the  equation  of  state  given  by  Eq.  (5-22)  is  used, 
Eq.  (H-3)  becomes, 

80D  +  28QMD  +  i   9Q(Du2)  +  |  (3oM)Du2 

+  4  [e-M/28r(Dur)  +  e"M/2(3  M  +  ^)(Du  ) 

r 

+  e"^2  1  8„(DuJ  ♦  e"M/2  (I,.Ht  ^^)(Du„) 


+  e"M/2  F-OT  V°V 

+  e_M/2  'FlInT  V)(DV]  "  °'        (H-H) 

In  order  to  find  suitable  expressions  that  might  replace 
the  correlation  functions  <id>  and  <T3Qd>  in  Eqs .  (9-^)  and 
(9-5),  the  following  procedure  is  used. 

If  Eq .  (H-^4)  is  multiplied  by  a  special  function 

h(T,8  t)  and  averaged  it  becomes 
o 

<h30D>  +  2<h(30M)D>  +  |  <h3Q(Du2)> 
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♦  I  «bDu»0,M)>  ♦  ^<e-M/2  F1±ni   »  (Du  )> 


-M/2   hDu,j, 


It  should  be  noted  that  in  Eq .  (H-5)  only  partial  deriva- 
tives with  respect  to  x°  and  4>  appear  explicitly.   Those 
terms  which  contain  partial  derivatives  with  respect  to  r 
and  0  are  buried  in  the  terms  containing  the  <j>  partial  deri- 
vatives.  To  see  how  this  can  be  done,  consider  the  follow- 
ing argument.   In  Eq .  (H-4),  there  are  three  types  of  terms 
involving  spatial  partial  derivative  with  respect  to  r,  0, 
and  $.      On  a  spatial  hypersurface  of  constant  x°  it  is 
reasonable  to  assume  that  the  average  of  a  quantity  associa- 
ted with  a  given  direction  in  space  ought  to  be  equal  to 
the  average  of  that  same  quantity  associated  with  a  different 
direction.   Since  the  r,  8,  and  <j>  terms  of  Eq.  (H-4)  have 
identical  roles,  (i.e.  transport  of  mass  in  r,  6,  <f>  directions 
as  Adler  [1]  points  out)   then  the  average  of  these  terms 
(multiplied  by  the  function  h  which  has  no  explicit  spatial 
dependence)  should  be  equal. 

Now,  D  can  be  decomposed  into  its  average  and  fluctuat- 
ing components.   That  is, 

D  =  p(x°)  +  d(x°,  r,  6,  <|>;0  (H-6) 

such  that  <D>  =  p"  . 

If  Eq.  (H-5)  is  expressed  in  terms  of  t  instead  of  M, 
it  becomes,  using  Eq.  (H-6) 
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90TT<h>   +    <h30d>   +   ~^-   p<h>   +   ^-<hd> 


+    2P<h 


V 


U+Tj 


R 


>    +    2<h 


V 


1+T 


d>    +   |   3Qp<hu2> 


+   |  P<h30(u2)>   +   |  <hu2O0d)>   +  |  <hd(8Qu2)> 


+  ^   P  ~  <hu2>   +  |   p<hu2 


+  I  g-  <hdu2>'+  |  <hdu2 


u,h 


V 


1  +  T 


V 


1  +  T 


+    l,<e-M/2   _X_   3    ?>    +    4-<e-M/2    _h 

r    sin    69  r    sin    6      cj>    <J> 


+  .<e-^2 


r  sin  6   4) 


3^d>  +  4<e 


-M/2    hd 


3,u  > 


r  sin  6   <j>  <{> 


+  ^<e"M/2 


hu 


* 


r  sin  6 


9^T 


+  H<. 


-M/2   hdu4) 


r  sin 


^1  +  T 
8  ,T 


1+T 


>  =  0  . 


(H-7) 


It  should  be  expected  that  weak  fluctuations  of  geometry 
do  not  create  pressure  gradients  which  accelerate  the  radia- 
tion fluid  to  local  velocities  near  the  speed  of  light. 
That  is, 

u2<<1 


whenever 


T*«l 


(H-8) 


This  makes  certain  moments  of  Eq.  (H-7)  small  compared  to 
others . 

Consider  now  two  cases.   Let  h  symbolically  stand  for 
two  special  functions,  t  and  9qt.   Then 
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<h>  =  0  . 


(H-9) 


Thus,  if  the  higher  order  moments  of  Eq .  (H-7)  are  neglected 
by  virtue  of  the  smallness  of  u2,this  equation  becomes 

<h9Qd>  +  — -  <hd>  +  2p<h30T>  +  |  p<h90u2> 


+  |  p<hu290x>  +  4p<e~M/2  5 — s-  3Au.> 

3       °  r  sin  8   <J>  <J> 


+  l\<t 


-M/2 


hu 


r  sin  6   (J) 


9,d>  +  kp<e 


-M/2    hu4> 


r  sin  9   4> 


3,t> 


=  0  . 


(H-10) 


As  pointed  out  by  Harrison  and  others  [9],  any  perturba- 
tion of  the  radiation  fluid  characterized  by  d,  u  ,  uQ,  u. 
will  propagate  with  the  sound  velocity  of  the  fluid.   For 
radiation  this  velocity  is  near  the  speed  of  light. 


V 


/dP  _  1 


sound    /  dD 


/3 


Hence,  such  an  accoustical  signal  would  propagate  along 
"near"  geodesic  lines.   That  is,  ds2  (accoustical  signal) 
is  near  enough  to  zero  that 


^  o     H/2,     M/2     .   ft  ,  . 
dx°  s  e   da  -  e    r  sin  0  d$ 


(H-ll) 


But  this  implies  that  the  three  moments  which  contain  the 

spatial  partial  derivatives  with  respect  to  <j>  in  Eq.  (H-10) 

-M/2 
ought  to  be  of  the  same  order  as  those  for  which  e      x 


1 


r  sin 


8cf  is  replaced  by  9  .   For  example, 


<e 


-M/2 


3,u  > 


r  sin  G   4>  $ 


<h3oV 
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Then,  Eq .  (H-10)  becomes  after  discarding  the  insignificant 
moments 

—  4R'    —     — 

<hd    d>  +  -^—  <hd>  +  2p<h8QT> 

If  the  first  two  terms  of  Eq .  (H-12)  are  consolidated 
into  one  term  f x<h8 od> (f or  an  appropriate  f 1  ) ,  it  becomes 

f  <hd    d>  +  ^p<e"M/2  5 — -  9  u  > 

10      p       r  sin  6   <J>  <J) 

=  -  2p<h8oT>  .  (H-13) 

If  the  characteristic  time  for  the  fluctuation  is  comparable 

1  -p 
to  the  characteristic  time  —  5-7-  of  the  average  geometry, 

then  the  second  term  is  equal  in  magnitude  to  the  first. 

In  this  case,  f . ,  which  is  a  function  of  x° ,  is  of  order 

one . 

Nov;  consider  the  remaining  equations  TJV;  v=  6  for  j 

equal  to  one,  two,  and  three.   After  some  rather  tedious 

calculations,  these  equations  become 

(Tlv;  v  =  0) 

8  [(D  +  £-)  e_M/2u  1  +  (f  3  M)[(D  +  — )  e_M/2u  ] 
0       c2         r      2   0         c2 

+  3  {e""M[Du*  +  L.(i  +  u*  -  u2)]} 

c2 

+  (|V,  +  1)  e-I;[Du^  +  ^  (1  +  u^  -«*)] 

r  c 

-  \   OrM    +   h    e"M[Du2    +  L_   (3   .    2u2)]  (H-l4) 

r  c2  cont. 
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p        ~M 


cont . 


c  p  -M 

+    (f  3    M   +   cote)[(D    +   — )  - —  u   u    1 

2      6  2  r  r   6 

c 


-M 

+  8,[(d  +  — )     e  .  n  u  uj 


<f> 


2      r    sine      r   <f> 


t-  p  -M 

+   |   OaM)[(D    +   — )    -£-r-fl 

2        $  2      r    sine      r   <j> 


u   ul    =   0  (H-14) 


(T2V;    v    =    0) 


P        e"M/2  5  P        e"M/2 

3»[(D   +  — )    ^-zr—  uj    +    (4  9   M)[(D   +  £-) 


r  6  2       o" 


uj 


r  0 


p  ~M 

+    3      [(D   +  — )    - —  u   uj 
r  2        r        re 

c 


5  4  P        e~M 

+   (t  3  M  +  ±)[(d  +  £-)   V-  u„un; 


+  ?  8e 


2    V         r 

f      -M 


r        r 


e_[Du'  +£_(1  +  u>  _ua}] 


2    r      6  r 


£_[Du>    +Z_(1   +   u>    .  u.}] 


I  <F  9eM) 


-M 

: [Du2    +  —   (3    -    2u2)] 

r  c2 


-M 


_   c^te    \    e_   [Dua    +   P_   (1   _   u2   _  u2)] 


>  c2 


-M 


K—r   9T(D   +   £-)    ^—  uQuJ 


r    sine      (j) 


r       e  <$>■ 


(H-15) 
cont . 
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+  (4 


5        1 


P    x    e 


-M 


cont . 


2  FXIKe  V)[(D  +  — )  V  uoV  =  °       ^-W) 


,      3  v 

(T      ;    v   =    0) 


-M/2 
3oC(D  +  77>  HlHS"  V   +   Cf  ?,M)[(D  + 


.    .     .5    ,   -xr^        P    ^    e~M/2 


2      r    sin6      c{) 


uj 


+    9r^D  +  ~)   Fl5nTurV 


,5        4        P     e"M 
+  (4  8  M  +  -)[(D  +  — )   e  .  Q  u  uft] 
2   r    r/LV      2   r  sine   r  6J 


1  P  ~M 

+  F  V(D  +  ^  fW  ueV 


— M 
4.    ^    1    .    M    ,    3    coteU/n        Pn      e  -, 

+  (2  F  V1  +  — F~ )[(D  +  ^  F^THe  ueV 

e"M  2  p  2 

FTiKe  %      F-5TE6  [%  +  ~2  (1  +  u*  "  u')] 


♦  (I      X 


2   r    sin6      <J> 


3.M) 


e~M  2         P  2        , 

■        ■    ft    [Du,    +  —   (1    -   u,    -u2  )1 
r   sine  4>  2  4> 


-M 


l> 


-  (lF-iinT  V}|  ^iIKe-[Du2  +IT  (3 


=  o   . 


(H-16) 


If  the  u   terms  are  discarded  and  if  the  equation  of 
state  is  used,  these  equations  become  simply 


ilRl  28oT 

VDur}    +  X%   +  TTF  (Dur' 


30(Du9)    +  -S-Du      ♦  -jyf   (Due) 


|  .-"V 


1      -M/2    1 

¥  e  F  8eD 


(H-1'7) 

cont . 
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llR  '  28  T 


1   -M/2    1 

¥  e 


9.D 


r  sin0   (j) 


(H-17) 


In  order  to  learn  something  about  the  behavior  of  the 

-M/2    h 
moment   <e  8(t)U(})>  wnlcn  appears  in  Eq .  (H-13), 

consider  the  average  of  the  third  of  Eq.  (H-17)  multiplied 


by   h. 


-^         ^  1  ,    -M/2        h 

p<h8nuA>    =   -  tj-  <e 


* 


r    sinG      cf> 


3,  d> 


(H-18) 


(To  obtain  Eq.  (H-18),  the  higher  order  moments  are 
neglec  ~ed . )  . 

Because  of  the  previous  arguments  concerning  the  rela- 
tionship of  spatial  and  temporal  partial  derivatives,  the 


following  equation  should  hold  for  some  f   of  order  one. 

2 


-  -M/2    h  „   .   ■*■ 

p<e     : — -  3,u,>  -  f  <hg  d> 

r  sine  $   §  2   °o 


(H-19) 


If  Eq.  (H-19)  is  substituted  into  Eq.  (H-13),  the  follow- 
ing equation  results. 


<h30d>  = 


-2 


fi  +  4f2 


p<h8.T>   . 


(H-20) 


If  the  changes  of  d  and  t  occur  over  the  same  time, 
it  must  be  true  that 


<hd>  =  f  p<hi> 

3 


(H-21) 


where  f   is  a  weak  function  of  x°  of  order  one  (unless  by 
chance  f,  +  4f   is  small). 

1         2 
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But  in  Eqs.  (9-4)  and  (9-5),  the  correlations  <Td> 
and  <O0x)d>  are  required.   Since  the  ratio  of  D  to  D  is 
(1  +  u2/3   )  (see  Eq.  (9-6))  which  is  very  close  to  one, 
it  can  be  assumed  that  for  some  functions  Xf  and  6f,  both 
weak  functions  of  x°  and  of  order  one  as  described  above, 
the  following  equations  hold. 

<id>  =  -  2Afp<i2> 

<(90T)d>=  -  3fp80<T2>.  (H-22) 

If  Xf  and  3f  are  expanded  in  a  Taylor  series  and  only 
the  constant  terms  are  retained  because  of  the  supposed 
weak  dependence  of  these  functions  on   x°,  the  following 
approximate  equations  result. 

<id>  =  -  2Xp<T2> 

<(3.T)d>=  _  3p90<t2>  •  (H-23) 

These  constants  are  then  to  be  viewed  as  arbitrary  parameters 
of  the  model.   Eqs.  (H-23)  are  therefore  consistent  with  the 
conservation  equations. 
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APPENDIX  I 

The  Dynamics  of  the  Energy  Momentum  Tensor 

for  a  Matter  Pilled  Universe 

The  dynamics  of  the  energy-momentum  distribution  is 
given  by 

T^v;  v  =  0.  (1-1) 

The  conservation  of  energy  equation  corresponding  to  y 
equal  to  zero  is  given  by  Eq .  (H-3)  with  u.  and  P  equal  to 

J 

zero.   (It  should  be  noted  that  u.  and  P  equal  to  zero  with 
D  unspecified  are  solutions  to  the  three  momentum  equations 
for  which  y  =  1,  2,  3.   See  Eqs.  (H-l4),  (H-15),  and  (H-l6).) 
That  is, 

80D  +  |(30M)D=  0    .  (1-2) 

3/2M 
If  this  equation  is  multiplied  by  e     and  integrated, 

the  following  equation  results. 

D  =  73S72      .  (I"3) 


With 


and 


D  =  p  +  d  (1-4) 


e3/2M  =  e3/2G  (1+T)3/2  f  (I_5) 

equation  (1-3)  becomes  for  small  t 

p  +  d  =_A_(1  _ix).  (1-6) 
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Hence , 

0     "           A 

P        e3G/2 

and 

3 

d    =   -   j  px 

(1-7) 


(1-8) 


That  is,  for  a  pressure  free  universe,  the  stochastic  process 
d  is  a  function  of  the  stochastic  prosess  t.   Hence, 

I   p<t2>  (1-9) 


and 


<xd>  =  -  5"  p<t^> 


<O0i)d>  =  -  f  p30<t2>   .  (1-10) 
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APPENDIX  J 
Solution  of  the  Growth  Equations 
for  a  Radiation  Filled  Universe 

2 


If  y  =  <t2>,  z  =  R-  <(^°T)     ,  x    R/RH 
Eqs.  (9-34)  and  (9-35)  are 


,  x  =  R/R  ,  the  growth 
ar 


x2  l!z+  x  dj  _  8Xy  =  2x^z  (J_1} 

dx 


2     dx 


and 


^+^z  =  ^  (J-2) 

dx   x      i,  dx  v    ' 

x 

+ 
where  R   is  a  constant  with  the  dimension  of  length.   All 

the  variables  (y,  z,  x)  in  Eqs.  (j-1)  and  (J-2)  are 

dimensionless . 

The  solution  of  these  two  coupled  differential  equations 

proceeds  as  follows.   If  Eq.  (J-2)  is  multiplied  by  x1*, 

the  left  hand  side  becomes  a  derivative  of  the  function  x'*z. 

This  equation  may  then  be  integrated  directly  to  give 

y  -  ^  x*z  +  fs  (J-3) 

where  f   is  an  arbitrary  constant  of  integration. 

o 

From  Eq.  (J-3),  the  first  and  second  derivatives  of  y 
with  respect  to  x  may  be  found  in  terms  of  z  and  x.   If 
these  are  substituted  into  Eq.  (J-1),  a  second  order  differ- 
ential equation  in  z  results. 
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dz  32pXf« 

—  +  9x  ~   +  8(2  -  A  -  B)z  =  i 

2       ax  i, 


d2z 


(J-4) 


dx 


Equation  (J-H)    can  be  transformed  into  a  second  order 
differential  equation  with  constant  coefficients  by  the 
change  of  variable 


Thus, 


x  =  e 


—  +  8  If-  +  8(2 

2     dt 


dt 


-  X  -  3)z  =  323Xf  e 

g 


-H 


(J-5) 


(J-6) 


The  associated  homogeneous  equation  of  Eq.  (j-6)  is 

—  +  8  ~  +  8(2  -  A  -  3)z  =  0  . 
dt2     dt 


Its  complementary  solution  z   is 


z   =  Aeit  +  BeJt 


(J-8) 


where    i  =  -  H    +  2/2(A+3) 
j  =  _  n  _  2/2(A+3) 

and  A  and  3  are  arbitrary  constants. 

The  complementary  solution  z   can  be  used  to  find  the 
general  solution  to  Eq .  (J-6).   This  is  done  by  assuming 
the  general  solution  is  given  by  Eq.  (J-8)  but  with  A  and 
B  functions  of  t. 

If  this  solution  form  is  substituted  into  Eq.  (J-6), 
a  differential  equation  in  A  and  b  results. 


A 


—  +-  8  ^   + 

dt2     dt 


(2  -  A  -  3) 


it 


(J-9) 
cont . 
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+  b  c  4i-  + 

dt2 


d_ 

dt 


+    8(2    -    A   -3)1    e 


Jt 


(J-9) 


r    dA   d(elt:)         dB   d(eJ,t)    , 
L    dt        dt  dt        dt        J 


+   |_[(e»)dA    +    eJtdB]+8[(elt)dA+eJtdB] 


32eAfge 


-kt 


(J-9) 


The  first  two  terms  are  zero  by  virtue  of  the  fact  that 
e   and  e    are  solutions  to  Eq.  (J-7).   In  the  equation, 


z(t)  =  A(t)eit  +  B(t)eJ,t 


(J-10) 


three  dependent  variables  of  t  are  present:  z,  A,  and  B. 
But  in  general,  three  variables  can  be  made  to  satisfy  three 
equations.   However,  so  far,  only  two  such  equations  have 
been  specified,  namely  Eq.  (J-4)  and  Eq.  (J-10).   Hence, 
a  third  equation  can  be  freely  chosen.   A  convenient  choice 

is 

it  dA  +   jt  dB  =  Q  (J_1;L) 

dt       dt 

so  that  Eq.  (J-9)  becomes 


d(elt)  dA   d(eJt)  dB  = 

dt   dt  +   dt   dt    ^$xlge 


-H 


(J-12) 


These  two  equations  can  be  written  in  matrix  form. 


Zu  =  F 


(J-13) 


where    Z  = 


it 


d(e   ) 
dt 


d(eJt) 
dt 
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u  = 


dA 
dt 

dB 
dt 


P  =  323Af  e 


By  Cramer's  Rule, 


-H 


f  _     N 


1 
v.     J 


dA 
dt 


det  Z 


det 


dB 
dt 


det  Z 


det 


323Af  e 
g 


it 


d(e   ) 
dt 


-H 


so  that 


dA 
dt 


32A3f  e 

_J5_ 


-(i+4)t 


J-i 


d(eJt) 
dt 


32ABf  e 


-4t 


(J-14) 


(J-15) 


dB 
dt 


32A3f  e 


-(J+4)t 


J-i 


(J-16) 


for     j  ^  i  (a  =  -3)  . 

Now,  Eqs .  (J-15)  and  (J-16)  may  be  integrated  to  give 


32x^fg    -(i+H)t   - 
(j-i) (l+4) 


(J-17) 


32A3f      (J  +  1|)t 
B  ~  "  (j-i)(j+4)  e       +  E 


(J-18) 


where  D  and  E  are  arbitrary  constants.   Hence,  by  Eq.  (J-10) 
and  by  the  identity  x  =  e  , 
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=  -  ^AFx""*1  +  Dx1  +  ExJ' 


(J-19) 


Sf 
where   F  =  t~- &■   and  is  a  new  arbitrary  constant.   By  Eq. 

(J-3)> 

y  =  F  +  gg  xi+ij  +  |g  xj+4  .  (J-20) 


If  these  solutions  for  y  and  z  are  substituted  into  Eqs. 
(J-l)  and  (J-2),  they  are  found  to  satisfy  these  equations. 
If  y  and  z  of  Eqs.  (J-19)  and  (J-20)  are  translated  back 


2 


2 


into  <x  >  and  <(90t)  >,  Eqs.  (J-19)  and  (J-20)  become 


<^>.=  Fl/3  +  Tf  Rl+,4+Tf  rJ+"  <j-21> 

2 

<(8oT)  >  =  -  i\\F1  /.R"^  +  I.  /.R1  +  Dn  7,RJ       (J-22) 
a_  1/3       1/3      1/3 


where    i  =  -  4  +  2/2(A+3) 

j  =  _  4  _  2/2(A+B) 

and  1-,/oi  ^i/v  an(^  ^i /^  are  arbitrary  constants.  (1-  .- 

D/R+   H  and  D  .  =   E/R+J   .) 
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APPENDIX  K 
Solution  of  the  Growth  Equations  for  a  Matter  Filled  Universe 

.3 

If  y  =  <t2>;  z  =  - —  <(3  t)2>;  x  =  —r   (for  some  R+) ,  the 

ad      o  r 

growth  Eqs.  (9-65)  and  (9-66)  are 

x2y"  +  |  xyf  -  3y  =  2x3z  (K-l) 

zf  +  £  z  -  !  i-  y«   .  (K-2) 

x3 

If  Eq.  (K-l)  is  differentiated  and  Eq.  (K-2)  is  substi- 
tuted into  it,  a  third  order  differential  equation  for  z 
results . 

2x3z,u    +    27x2zn     +   84xz'    +    42z    =    0    .  (K-3) 

If  a  change  of  variables  is  made  such  that 

x  =  et  (K-4) 

then  Eq.  (K-3)  becomes 

2z  tt  +  21z  t  +  6lzt  +  42z  =  0  .  (K-5) 

The  general  solution  of  this  equation  is  simply, 

z  =  Ge^   +  lde~lt/2   +   Je"6t   .  (K-6) 

In  terms  of  x,  Eq .  (K-6)  becomes 

z  =  Gx"1  +  Hx~7/2  +  Jx~6  .  (K-7) 

If  this  equation  is  substituted  into  Eq.  (K-2),  an  equation 
for  y'  results.   If  this  equation  is  integrated,  the 
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following  expression  for  y  appears. 

y  =  G  x2  -  |S  x"1/2  +  |  jx-3  +  k  .  (K-8) 

But  Eqs.  (K-7)  and  (K-8)  are  general  solutions  to  the 
set  of  the  two  coupled  Eqs.  (K-3)  and  (K-2).   Four  integra- 
tions are  required  to  solve  this  system  which  yields  four 
arbitrary  constants  (G,  H,  J,  and  K  above).   However,  the 
given  system  of  equations  (Eq.  (K-l)  and  (K-2))  require  only 
three  integrations.   Because  Eq.  (K-3)  is  derived  from  these 
equations  by  differentiation,  Eqs.  (K-7)  and  (K-8)  must  be 
the  general  solution  to  the  given  equations  with  one  of  the 
appropriate  arbitrary  constants  equal  to  zero.   By  inspection, 
the  equation 

y  =  K  and  z  =  0 

although  a  solution  to  Eqs.  (K-3)  and  (K-2),  is  not  a  solu- 
tion to  the  original  equations.   Hence,  the  general  solution 
to  Eqs.  (K-l)  and  (K-2)  must  be 

z  =  Gx_1  +  Hx~7/2  +  Jx"6  (K-9) 

77  2    2H   -1/2    4  T  -3  (v    ins 

y  =  Gx   -  ~—  x      +q-Jx.  tK-10) 

Hence , 

2 

<^oT)  >  _i       -7/2       -6 

=  I  R    +  FnR  u      +   DnR  u  (K-ll) 


and 


a,       o       o         o 
d 


<T2>  =  I0IT  -  —  R  1/d    +  ^  PQR  5     .  (K-12) 
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APPENDIX  L 
Classification  of  Solutions  for  the  Radiation  Filled  Universe 

2 

With  the  exception  of  the  integration  of  R'  ,  Eqs .  (E-20) 
and  (E-21)  of  Appendix  E  give  the  solutions  to  Eqs.  (E-l) 
and  (E-2).   These  solutions  are,  with  k  neglected  as  usual 
(see  Eqs.  (9-4l)  to  (9-^5)), 


where 


Kp  =  _  3^.  f  RmA 
R,    A 


R'2  =arR~2  +  bARmA+2 


f. 


1/3 


D 


1/3 


C(Il/3) 
12 


C(Dl/3) 

12 


10X  +  43  -  1  + 


10X  +  43  -  1  - 


23 


/2(X+3) 


23 


/2(X+3) 


1/3 


C(Il/3) 
~3^ 


2X  +  /2(X+g) 
/2(A+3) 


'Dl/3)  f  2X  -  /2(X+3) 


D 


1/3 


!£ 


/2(X+3) 


m. 


=  i  =  -  4  +  2/2(X+3) 


1/3 


mn     =  j  =  _  l|  _  2/2(X  +  3)   . 
Ul/3, 


(L-D 
(L-2) 

(L-3) 
(L-4) 

(L-5) 

(L-6) 
(L-7) 
(L-8) 


The  A  subscript  designates  the  type  of  fluctuation.  (I 


i/: 


type  corresponds  to  a  situation  where  the  nonunif ormity 
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increases  with  increasing  R  and  D  ,_  -  type  corresponds  to 
a  situation  where  the  nonunif ormity  decreases  with  inceasing 
R.) 

There  are  four  kinds  of  solutions  based  upon  the  manner 
in  which  Kp  and  R  are  displaced  from  the  fluctuation  free 
solutions  Kp*  and  R* .   That  is,  f.  and  b   can  in  principle 
be  each  positive  or  negative  depending  on  the  choice  of  A 
and  3. 

In  order  to  determine  how  a  given  (A, 8)  displace  Kp  and 
R  from  Kp*  and  R* ,  it  is  necessary  to  locate  those  values 
(A, 8)  for  which  f  and  b  are  equal  to  zero.   Then  a  set  (A, 8) 
which  lies  on  a  particular  side  of  the  two  curves  of  zero 
f.  and  zero  b.  determines  the  sign  of  f.  and  b.  and  hence 
the  direction  of  displacement  of  kp  and  R  from  Kp*  and  R* . 

In  order  to  find  the  line  of  zero  f.  in  A-8  space,  con- 
sider first  the  following  definition. 

n  =  /2(A+8)   .  (L-9) 

The  set  (A, 3)  for  which  f.  equals  zero  is  given  by  the 
equation 

10A  +  48  -  1  =  +   23  (L-10) 

/2(A+8) 

where  the  plus  sign  is  used  if  A  represents  D-./o  and  the 
minus  sign  if  A  represents  I-,/?-   This  equation  as  it  stands 
is  rather  difficult  to  plot  directly.   But  the  plotting  is 
simplified  if  Eq.  (L-10)  is  separated  into  two  parametric 
equations  with  parameter  n. 
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The  two  parametric  equations  3(n)  and  X(n)  are  found  to 


be 


R  =  - 
3   2 


X  =  -  - 
A     2 


5n2  -  1 
3n  +  1 

2n2  + 


(L-ll) 

(L-12) 
type 


3n  +  1 

where  n  >  0.  The  top  sign  corresponds  to  the  D  /0 
fluctuation  and  the  bottom  sign  to  the  1-,/n  -  type 
fluctuation. 

By  varying  n,  the  equation  3(A)  can  be  traced  out.   This 
equation  is  plotted  in  Figs.  3  and  4  for  both  I-.  •,  and  D..  ,_ 
type  fluctuations  respectively. 

The  equation  of  X  vs.  3  for  b   =  0  is  much  simpler  to 
find.   With  straightforward  algebra,  this  equation  is 


3  =  X(2X  -  1) 


(L-13) 


for  both  the  I,  ,~   and  D,  ,~  -  type  fluctuations.   This  equa- 
tion is  also  plotted  in  Figs.  3  and  4.   It  should  be  noted 
that  Eq.  (L-13)  is  quadratic  in  X.   Hence,  according  to  this 
equation,  there  should  be  two  values  of  X  corresponding  to 
a  given  value  of  3.   However,  only  one  value  satisfies  the 
requirement  that  b.  equal  zero.   That  there  are  two  values  of 
X  comes  from  the  fact  that  Eq .  (L-13)  is  derived  by  a  squar- 
ing operation  which  thereby  introduces  an  extra  root. 
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APPENDIX  M 
Index  and  Location  of  Symbols  Used  in  Dissertation 

The  following  list  contains  the  symbols  used  in  this 
work.   Accompanying  each  symbol  is  the  location  where  the 

definition  or  explanation  of  the  symbol  may  be  found. 

1.  L*,  G,  c  I 

2.  ds,  gyv,  xy,  x°,  r,  6,  d,  da2,  RQ ,  IV-A 

G    M 

RM,  k,  e  ,  e  ,  M  IV-A 

3-   f   ,  fM»  T  IV~B 
y  v    M 

**.   [yv,3J58  ,  ra  ,  gaB,  R   ,  RX  ,  R, 
'   5  y    yv        yv    v 

GX  ,  k  V-A 

5.  D,  P,  u,  D,  P,  a  V-F 

6.  R  VI-A 

7.  p,  d  IX-A,  B, 

8.  M,  V,  AD,  AM,  AV,  g,  V,  A,  3  IX-A-2 

9.  F1/3,  I1/3,  D1/3,  i,  j  IX-A-3 

10.  C(I1/3),  C(D1/3),  C(A)  IX-A-H-a 

11.  p*,  R*  IX-A-4-b 

12.  F  ,  I  ,  D  IX-B-3 

13.  B(Iq),  B(Dq)  IX-B-4 

14.  6p,  D,  I,  Tr  TB,  <i2>i5  <t2>0,  R  X-A 

13^ 


15 
16 
17 
18 
19 

20 
21 


6T,  T,  HR,  L,  6p  .  p  .  z,    R  .  R    .  R 


lR 


v,  V,  v,  P,  p 


o  '  r  o 


o J   max '   e 


F(6;xy),  f(D;xy)s  h[D(xy)] 
TyV   TA    Uy   U    y 


Pv,  Py,  av,  V,  in,  k,  (»),  C,  B,  a,  ap , 


a,,  f,  b 


p,  a,  h,  fj,  f2,  f3,  x   e 


f»  Hf 


y,  x,  z,  R  ,  t 


22.  t,  A,  B,  Z,  u,  F,  fg,  D,  E,  P 

23.  G,  H,  J,  K 


X-B 

A 
B 
D 

E 
H 

J,  K 
J 
K 


24.   m  ,  f     ,  f     ,  b     ,  b     ,  m 

A     l1/3  U1/3     I1/3    U1/3     l1/3 

m     ,  f  ,  b   n 

u1/3        A  a 
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Figure  1.   Classification  Designations 
for  the  Increasing  Type  Fluctuations. 
The  solid  lines  are  the  fluctuation  free  solutions  for  R 
and  p.   The  broken  lines  depict  the  possible  deviations  of 
R  and  p  from  these  unperturbed  solutions.   The  letters  K, 
L,  M,  and  N  label  the  four  possible  combinations  of  perturbed 
solutions  for  R  and  p  from  their  unperturbed  values. 
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Figure  2.   Classification  Designations 
for  the  Decreasing  Type  Fluctuations. 
The  solid  lines  are  the  fluctuation  free  solutions  for  R 
and  p.   The  broken  lines  depict  the  possible  deviations  of 
R  and  p  from  these  unperturbed  solutions.   The  letters  K,  L, 
M,  and  N  label  the  four  possible  combinations  of  perturbed 
solutions  for  R  and  p  from  their  unperturbed  values. 
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Forbidden 
Region 
(Exponent  i  is 
imaginary) 


^ 


~H/2        -I 


1/2 


1/2 


1 1/5 


Figure  3.   Solution  Types  in  A-$  Space  for  the  Increasing 
Type  Fluctuations  for  the  Radiation  Filled  Universe 
The  x-3  plane  is  sectioned  into  four  regions.   The  perturbed 
solutions  for  each  region  for  R  and  p  are  depicted  as  broken 
lines.   The  fluctuation  free  solutions  are  represented  as 
solid  lines. 
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Forbidden 
Region 

(Exponentj  is 

imaginary) 


-1/2 


Figure  4.   Solution  Types  in  X-$  Space  for  the  Decreasing 
Type  Fluctuations  for  the  Radiation  Filled  Universe. 
The  A-3  plane  is  sectioned  into  three  regions.   The  perturbed 
solutions  for  each  region  for  R  and  p  are  depicted  as  broken 
lines.   The  fluctuation  free  solutions  are  represented  as 
solid  lines. 
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Figure  5.   Lines  of  Constant  i  and  j 


mo 
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